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SUMl^IARY 

The  present  report  sxxmmarlzes  the  results  of 
a  study  which  the  New  York  University  Group  of  the 
Applied  Mathematics  Panel  has  undertaken  upon  the 
request  of  the  Bureau  of  Aeronautics,  Navy  Department 
(Project  No.  llA.j;167),  and  which  was  carried  out  under 
the  responsibility  of  Professor  K.O.Priedrichs  with 
the  assistance  of  Dr.  Chas.  R.  DePrima  and  other  mem- 
bers of  the  group.  Our  group  is  greatly  indebted  to 
Mr.  E.  S.  Roberts  of  the  American  Cyanamid  Corr5)any 
for  his  stimulating  advice. 

The  original  request  was  for  assistance  in  the 
design  of  unconventional  nozzles  for  rocket  motors, 
with  a  view  toward  attaining  shortness  while  preserving 
high  efficiency.  Since  the  opening  angle  of  such 
nozzles  is  necessarily  wide  and  the  expansion  rapid, 
it  was  necessary  to  improve  the  classical  (hydraulic) 
theory  of  the  de  Laval  nozzle  and  to  widen  the  scope 
of  the  investigation  further  by  consideration  of  gas 
dynamical  phenomena  more  generally.  Therefore,  the 
present  report  is  not  concerned  solely  with  the  original 
problem  of  exhaust  nozzles,  but  it  contains  methods 
and  results  of  potential  interest  for  other  problems j  in 
particular,  attention  is  given  to  the  design  of  "kinetic 
compressors",  i.e.,  of  nozzles  serving,  not  as  means 
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for  ejecting  a  supersonic  Jet,  but  for  receiving  a 
parallel  flow  of  gas  at  high  supersonic  speed  in  order 
to  compress  and  arrest  the  incoming  gas» 

More  specifically,  the  contents  of  the  present 
report  can  be  summarized  as  follows; 

Part  I     contains   a  three-dimensional   treatment 
of  flow  through  nozzles  as  a  refinement  of  the  customary 
one -dimensional  theory,  under  the  following  ass\iraptions: 

1)  The  fluid  is  ideal  and  homogeneous,  and  the  flow 
isentropic,  steady  and  irrotational.  2)  Viscosity  and 
heat  conduction  are  ignored.   3)  It  is  further  assumed 
that  shocks  and  jet  detachment  do  not  occur  in  the 
nozzle;  however,  conditions  for  the  absence  of  these 
phenomena  are  analyzed* 

A  simple  formula  for  the  thrust  produced  by  such 
a  flow  through  a  nozzle  is  derived  and  applied  to  various 
types  of  contours;  particular  emphasis  is  given  to  widely 
divergent  short  nozzles.   Among  families  of  nozzles  yield- 
ing the  same  thrust,  the  shortest  nozzle  is  determined. 
Tables  I  to  III  at  the  end  of  Part  I  (p.  35)  show  numeri- 
cal results  in  a  condensed  form. 

Part  II  concerns,  first,  the  construction  of  "perfect" 
exhaust  nozzles  which  yield  maximum  thrust  and  a  parallel 
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exhaust  flow  for  a  prescribed  expansion  ratio.  Secondly, 
reverse  nozzle  flow  or  "compressor"  flow  Is  discussed. 
"Perfection",  which  is  of  little  Importance  for  exhaust 
flow,  is  here  essential. 

The  stability  of  the  conpressor  flow  is  then 
discussed;  it  is  shown  that  shocks  must  be  admitted  to 
avoid  instability  of  the  phenomena  under  slight  varia- 
tions of  velocities  and  pressures. 

Furthermore,  conditions  are  analyzed  for  the  possi- 
bility of  a  supersonic  stream  of  air  entering  a  coitrpressor 
without  interference  of  shocks.  Limitations  of  the  rim 
angle  are  given  for  a  compressor  carried  in  the  nose  part 
of  a  projectile. 

Finally,  some  remarks  on  the  drag  against  a  compressor 

carrying  projectile  are  added. 

The  Appendix  supplies  mathematical  detail  for  Part  I. 
Furthermore,  it  contains  an  exposition  of  the  mathematical 
procedure  for  the  construction  of  perfect  nozzles,  the 
computational  schemes  developed  in  the  Appendix  can  easily 
be  applied  to  other  problems  of  a  similar  character. 

R.  Courant 
Contractor»s  Technical  Representative 
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THSORETICAL  STUDIES  ON  THE  FLOW  THROUGH  NOZZLES 
AND  RELATED  PROBLEMS 

Nozzles  of  customary  design,  ending  In  a 
cone  of  small  angle  of  opening,  have  proved  to  be 
very  satisfactory  for  the  purpose  of  rocket  propul- 
sion.  Thrusts  between  90  percent  and  100  percent  of 
the  value  predicted  by  gas  dynamical  theory  were 
foiond  in  a  great  number  of  tests.  Also,  when  the 
nozzle  was  somewhat  shortened  or  lengthened  or  if  the 
angle  of  opening  was  varied  by  several  degrees,  the 
resulting  thrust  remained  quite  satisfactory.  This 
means  that  not  much  increase  in  efficiency  can  be 
expected  from  refinements  in  nozzle  construction* 

There  are  situations,  however,  calling  for 
deviation  from  conventional  nozzle  design.   In  partic- 
ular it  may  be  desired  to  make  the  nozzle  rather  short 
and  accordingly  the  angle  of  opening  rather  wide.  For 
such  nozzles  the  usual  "one -dimensional"  or '^hydraulic" 
thebry  is  no  loiter  adequate  and  a  more  refined  treat- 
ment is  necessary  to  obtain  theoretical  information. 
Such  a  refinement  of  the  hydraulic  treatment  is  developed 
in  the  first  part  of  this  memorandum.   Simple  formulas 
are  derived  for  pressure  distribution  and  thrust.   These 
formulas  are  more  accurate  than  those  derived  from  the 
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hydraulic  theory  and  are  therefore  applicable  to 
shorter  nozzles  with  somewhat  wider  angles  of  opening. 

According  to  the  hydraulic  theory,  the  thrust 
produced  by  the  flow  through  a  nozzle  depends  only  on 
the  expansion  ratio,  i.e.,  the  ratio  of  exit  to  throat 
cross-section  area;  thus  short  nozzles  would  yield  the 
same  thrust  as  long  ones  if  their  expansion  ratio  is 
the  same.  One  cannot  expect  this  conclusion  to  be 
correct,  and  indeed  the  refined  theory  shov/s  that  the 
optimal  thrust  can  be  obtained  only  by  very  long  nozzles. 
Nevertheless,  it  remains  true  that  one  can  obtain  a  few 
percent  less  than  the  optimal  thrust  by  much  shorter 
nozzles. 

To  obtain  the  optimal  thrust  requires  a  partic- 
ular nozzle  design  which  we  shall  refer  to  as  "perfect". 
The  end  section  of  such  a  perfect  nozzle  bends  gently 
inward  so  as  to  produce  a  cor^jletely  axial  thrust 
flow  with  constant  velocity.  While  it  is  of  no  practi- 
cal importance  to  make  the  nozzles  perfect  if  their 
purpose  is  to  create  thrust,  it  is  necessary  to  design 
perfect  nozzles  for  other  problems  in  which  straight- 
ness  and  constancy  of  the  exhaust  flow  are  more  relevant 
than  shortness.   Apparently  this  is  the  case  for  super- 
sonic wind  tunnels  and  for  superaonie  diffusers  or 
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compressors.   A  compressor  is  a  perfect  nozzle  to  be 
operated  in  reverse  direction.   The  air  is  to  enter 
the  compressor  with  parallel  and  constant  supersonic 
velocity  and  is  to  be  compressed  past  the  throat,  poss- 
ibly so  as  to  reduce  the  velocity  to  zero* 

Of  particular  interest  are  supersonic  compressors 
which  are  carried  by  a  projectile-shaped  body  opposed  to 
a  uniform  stream  of  air.  The  problem  arises  of  deter- 
mining under  which  circumstances  the  air  would  enter  the 
compressor  with  constant  and  parallel  velocity  and  with- 
out interference  of  shocks. 

A  further  problem  is  that  of  the  stability  of  the 
flow  through  a  compressor.   It  will  be  shown  that  the 
"perfect"  shockless  flow  is  unstable  while  the  flow 
appears  to  be  stable  when  certain  shock  fronts  in  the 
interior  of  the  compressor  are  admitted. 

The  problems  of  the  perfect  nozzle  and  of  the 
supersonic  compressor  are  discussed  in  the  second  part 
of  this  memorandum. 

Detailed  mathematical  developments  are  placed  in 
the  Appendix, 
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PART  I. 
ISENTROPIC  FLOW  THROUGH  N0ZZL3S 
1,  Basic  Relations  and  Hydraulic  Treatment. 

Ihe  nozzle  to  be  Investigated  Is  a  "de  Laval" 
nozzle;  i.e..  It  is  of  the  "converging-diverging" 
type  and  has  axial  symmetry.   The  gas  is  to  leave  the 
chamber  (c)  with  velocity  zero,  to  become  supersonic 
while  passing  through  the  throat  (t),  and  to  enter  the 
atmosphere  when  leaving  the  nozzle  exit  (e).  Of.  Fig.  1« 
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We  assume  the  gas  to  be  ideal,  homogeneous^  "■'and  isentropic 
and  the  flow  to  be  steady  and  irrotatlonal.   Viscosity  and 
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•;:-  For  exhaust  gases  of  a  rocket  motor  this  assumption 
implies  that  the  fuel  has  burnt  completely  in  the  chamber; 
if  the  gas  is  not  homogeneous,  one  may  employ  an  average 
adiabatic  exponent  to  some  approximation. 
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heat  conduction  are  ignored.   This  appears  to  be  justi- 
fied for  the  exhaust  flow  (Cf.  T.  Stanton  [  1]  )^''""^ except 
that  viscosity  will  play  a  role  in  determining  where  Jet 
detachment  and  shocks  occur.  The  formulas  to  be  given 
are  valid  as  long  as  the  flow  remains  free  of  shocks 
and  as  long  as  the  jet  remains  attached  to  the  nozzle 
wall* 

The  isentropic  character  of  the- gas  flow  is  ex- 
pressed by  the  relation 


(1.01) 


p  ^   =  const. 


p  being  the  pressure,  x>  the  density  (mass  per  unit  volume) 
and  f  the  adiabatic  exponent.  The  sound  velocity 


(1.02) 


0  =  (  f  p  ^  -1)V2 


is  connected  with  the  flow  speed  q  through  Bernoulli's 
law  which  we  write  in  the  form 


(1.03) 


l-^)c  +^q  =c^=q^ 


*  Numbers  in  brackets  refer  to  the  bibliography  at  the 
end  of  this  memorandum. 
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The  significance  of  the  "critical"  speed  q^.  =  c,^ 
Is  that  the  f lov?  speed  Is  critical  when  it  coincides 
with  the  local  sound  speed*   Subsonic  and  supersonic 
flow  can  then  be  distinguished  by  q  <  q^  and  q  >  q^ 

since,  as  Is  easily  seen,  q  >  q^  implies  q  >  c  • 

With  the  aid  of  Bernoulli's  law  one  can  express 
the  quantities   c,  p,/o  in  terms  of  q/q^.  ; 
using  the  abbreviation 


V  = 


^  -  1 


(1.04) 


(1.05) 


(1.06) 


V 


1  -  /x2(q/q,,)^ 


»  (°.^.  =  <u) 


1  - 


/ 


2V  + 

2 

1   -    /i^(q/q^)^ 

e 

1   - 

>M^(q/q*)^ 

• 

e* 

n 

_    ..2 

V  +  1 
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The  "one-dir.ensional"  or  hydraulic  theory  of 
0»  Reynolds  (1885)  assumes,  as  a  first  approximation, 
that  the  flow  speed  q  and  hence  also  c,p,^  are 
constant  over  cross -sections  perpendicular  to  the 
nozzle  axis.   The  mass  flux  per  unit  time  across  a 
cross-section  of  area  A  is  then  given  by  A  /o  q;  hence 
this  quantity  is  a  constant,  or 

A      <°-«-^* 


/oq 


A.,  being  the  area  of  the  cross-section  at  which  the 
critical  speed  is  assumed.   One  easily  derives  the 
relation 


(1.07) 


with 


1,S^  =  (1  -  m2)  ia  , 
pq         q 

M^  -  1  =  (q/c)^  -  1  = 


2 


4 


2  2 
->Aq 


(m  =  q/c     being   the  Mach  number)  from  which  follows   the 
hydraulic  approximation. 


(1.08) 


M  =    (M^    -   1)  ifl      . 
A  q 
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This  relation  shows  that  a  flow  with  Increasing 
speed  q  is  possible  for  q  <  q^,.  or  q  <  c  only 
when  A  decreases,  and  for  q  >  q^,.  or  q  >  c  only 
when  A  increases.   In  other  words,  in  the  hydraulic 
approximation,  flow  with  increasing  speed  is  possible 
only  when  the  critical  speed  is  just  reached  at  the 
throat. 

2.  Refined  Treatment 

In  order  to  refine  the  hydraulic  treatment  just 
explained  it  is  necessary  to  employ  the  basic  partial 
differential  equations  for  an  irrotational  isentropic 
flow.  We  introduce  coordinates;  x,  along  the  axis  in 
the  exhaust  direction,  and  y,  the  distance  from  this 
axis;  further  we  introduce  the  angle  6    of  the  flow 
direction  with  respect  to  the  x-axis,  the  potential 
function  ^    and  the  stream  function  i/r   • 

The  stream  function  Y     ^^   to  be  so  normed  that 

2 
tr  y  is  the  rate  of  mass  flow  per  unit  time  carried  by 

the  stream  tube  yr   <  const.   The  basic  equations  for  the 

four  dependent  variables  ^,  ^^ ,  q,  and  6  can  then 

be  written  in  the  form 


(2,01)       ^  =  q  cose  ,^  =  q  sinG 
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(2.02)  -W^^  =  -  /oqy  ^^^^    »   '^^^  =  ^^^  °**^®  • 

The  density  p     Is  to  be  expressed  in  terms  of  q 
by  means  of  relations  (1.06)  and  (1.04). 

It  would  be  natural  to  prescribe  nozzle  contour 
and  state  In  the  chamber  and  to  ask  for  the  resulting 
flow.  Mathematically  this  would  mean  solving  a  boundary 
value  problem  for  equations  (2.01)  and  (2.02).  To  sim- 
plify the  task  we  reverse  the  procedure:  We  first  pre- 
scribe the  velocity  distribution  along  the  axis.  - 

(2.03)        q  =  %M   f or  y  =  0 

and  then  determine  possible  nozzle  contours  from  the  re- 
sulting stream  surfaces.  To  this  end  we  Introduce  the 
stream  functions   ^  and  ^f/    as  Independent  variables  in- 
stead of  X  and  y,  consider  the  quantities  tl,   y,  © , 
and  q  as  dependent  variables,  and  develop  them  with 
respect  to  powers  of  yi*^  It  is  sufficient  here  to  report 


*  One  may  also  Introduce  q  and  V  as  independent  vari- 
ables: this  might  seem  to  offer  advantages  since  the  non- 
linearity  would  then  refer  only  to  the  independent  variable 
a.  Nevertheless  it  was  found  that  the  present  scheme  is 
better,  at  least  for  the  expansions  up  to  the  second  order 
as  given  here. 
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only  the  results.   Details  will  be  carried  out  in  the 
Appendix,  Section  12. 

To  fonaulate  the  results  we  introduce  the  impor- 
tant dimensionless  quantity  h  by 


(2.04) 


h  =  ^f^^q^TJq. 


the  significance  of  which  is  that  h  =  A/A^,_  in  hydrau- 
lic approximation.   Since  by  (1.06)  the  density  a     is 
a  given  function  of  q  the  same  is  true  for  hj  we  have 
h  =  1  for  q  =  q.,..   Frora  (1.07)  we  obtain 


(2.05) 


2  dh/h  =  (M^^  -  1)  dq/q  . 


On  the   axis,   where  we  had  assximed     q  =  q    (x)      the    quantity 
h     becomes   also  a   function  of     x 


(2.06)   h  =  h^{x)  =  (q,,yq) 


1/2 


1  -  }-^Wi.,f 


1  - 


^' 


-2 


[Cf.  (2.04),  (1.06),  (1.04)]  which  is  given  inasmuch  as 
q.QM    is  given.   Also  M,  defined  by  (1.08),  becomes  a 
given  function  M^(x)  along  the  axis  when  q  =  q  (x)  is 
inserted. 
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We  place   the   origin     x  =  0,   y  =  0,    at  the 
throat;   more  precisely,  we  place   it  such  that 


(2.07) 


^o^O)   =  q, 


or  what  is  the  same  thing,  such  that  l^'(O)  =  0.  We 
then  have 


(2.08) 


h^(0)  =  1,  and  M  (0)  =  1  . 


Instead  of  0   and  '\^  ,  two  other  parameters,  | 
and  07   ,  are   introduced  by  the  relations 


(2.09) 


o 


(2.10) 


^=  /e^v 


Clearly,   ^   and  Oj  have  the  dimension  of  a  length  and 
\      reduces  to  x  on  the  axis  V  =  0,(ji5  =  0  for  x  =  y  =  0 
being  assumed)* 

In  terms  of  the  quantities  just  introduced,  the 
expansions  of  x,  y,  q,  6  ,  h  and  p  with  respect  to  powers 
of  >?  (instead  of  y   )  are.  If  h^(  ^  ),  h^(  ^  ),  hj(  \  ), 
Mq(  I  ),  and  Pq(  |  )  are  written  singly  h^,h^,hJJ,MQ,  and 
p^  respectively: 


JllUi 


(2.11)  ::  =  ^  -i  Vi'^^ 

(2.12)    y  =  h^77  1 1  +  1  [  (m2-  1)  h^h;;  -  (h.)^],^^! 

(2.13)  q  =  q^  jl+^h^hj  07  ^j 

(2.14)  0  =  h»>j 

(2.15)  h  =  h^  I  1  +  ^  (M^  -  1)  n^h;;  7,2| 

(2.16)  p  =  p^  {i-^M^h^h;;  7i^\ 

(For  a  detailed  derivation  see  Section  12). 
It  is  to  be  noted  that  the  siirfaces   ^  =  const. 
and  y}    ■  const,  are  the  potential  and  stream  s\irfaces 
respectively.   Thus  the  two  equations  (2.11)  and  (2.12) 
yield  parametric  representations  for  potential  and 
stream  surfaces.   If  q  (  ^  )  is  chosen,  and  then  h  (  ^  ) 
and  Mq(  I  )  are  determined  from  (2.05)  and  (1.08),  these 
stream  and  potential  surfaces  are  easily  drawn.  Each 


stream  surface  may  serve  as  nozzle  contour. 

2 

If  the  terms  involving  07  were  neglected  one 

would  obtain  x  =  ^  '^"^o^^^''?'  ory  =  h(x)07  .  This 
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relation  Is  identical  with  the  result  of  the  hydraulic 

theoryj  indeed,  since  x  =  0   corresponds  to  the  throat 

2   2 

and  h  (0)  =  1  we  would  have  A/A^^  =  (y/y..^)  =  ^q^^^' 

It  Is  thus  clear  that  the  formulas  above  represent  a 
refinement  of  the  hydraulic  theory. 

The  relation  (2.13)  for  q  yields  an  improvement 
over  the  hydraulic  assumption  that  the  speed  is  constant 
on  the  potential  surfaces.   To  form  an  idea  about  the 
magnitude  of  the  deviation  from  constancy  we  may  intro- 
duce the  radius  of  curvature  R  of  the  stream  contour 
*>)  =  const.   R   is  obtained  approximately  from 
y  =  h^(^  )<>]  ,  X  =  ^  ,  as 

R  =  l/hj'oi. 

Therefore,  relation  (2.13)  can  be  written  in  the  approximate 
form 


(2.17) 


q  =  q. 


,(i*H) 


This  formula  shows  that  the  deviation  of  the  speed  from 
that  at  the  axis  depends  on  the  ratio  of  the  width  of  the 
nozzle  to  the  radius  of  cxirvature  of  the  nozzle  contour. 

Unless  the  end  section  of  the  nozzle  curves  inward 
as  would  be  the  case  for  h^(  ^  )  <  0,  the  speed  Increases 


If'  ii  f"  ^ 


3 
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and  thus,  on  moving  away  from  the  axis,  the 
pressure  decreases  along  the  potential  sturfaces* 
This  is  equivalent  to  saying  that  the  lines  of 
constant  speed  and  hence  the  lines  of  constant 
pressure  bend  backwards  as  shov/n  in  Table  I, 
Our  formulas  thus  indicate  that  the  pressure  along 
the  nozzle  wall  assumes  ^Blues  that  are  less  than 
those  calculated  from  the  hydraulic  theory*   This 
behavior  could  be  expected  since  a  widely  divergent 
opening  will  give  the  gas  an  opportunity  for  quick 
expansion. 

3.  Critical  Remarks 

Before  deriving  formulas  for  the  thrust  and 
before  considering  exanples,  we  shall  discuss  several 
possible  objections  to  the  preceding  method. 

On  mathematical  grounds  it  may  be  considered 
improper  to  prescribe  a  quantity  such  as  q  (x)  along 
the  axis,  as  we  did,  whereas  the  proper  procedure 
would  have  been  to  prescribe  the  nozzle  contovir  and 
the  state  in  the  chamber,  in  agreement  with  what  is 
prescribed  in  reality.  A  way  of  prescribing  data 
mathematically  is  said  to  be  proper  if  for  such  data 
there  exists  a  unique  solution,  which  depends  continuously 


^^mmm-kmil^ 
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on  these  da ta»  'Prescribing  the  quantity  qQ(^)  on 
the  axis  is  indeed  not  proper  in  this  sense.   It  is 
true,  if  q  (x)   is  an  analytic  function,  a  unique 
solution  of  the  differential  equations  exists  in  the 
neighborhood  of  the  axis  but  one   does  not  know  how  far 
this  solution  extends  without  developing  singularities. 
It  is  also  known  that  if  <Iq(^)      i^  analytic  the  ex- 
pansion employed  converges  in  a  neighborhood  of  the 
axis  but  one  does  not  know  how  large  the  range  of  con- 
vergence is.  Further  it  is  known  that  slight  changes 
in  the  character  of  the  function  q^C^)  ^^7  cause 
considerable  changes  in  the  solutions  not  only  in  the 
subsonic  range  but  in  the  supersonic  range  as  well. 
This  will,  however,  not  occur  if  the  change  in  q^W 
is  "smooth"  and  analytic.  Our  method,  therefore,  seems 
to  be  justified  as  long  as  the  function  (IqM      is 
chosen  smooth  and  analytic.   The  same  character  is  then 
imparted  to  the  streamlines  in  the  neighborhood  of  the 
axis.   Thus  we  can  say:   Our  method  for  solving  the 
differential  equations  of  irrotational  isentropic  flow 
(2.01),  (2.02),  is  applicable  if  the  desired  nozzle 


«-  For  this  notion  see  [2]  p.  171-179.  That  prescribing 
nozzle  contours  and  state  in  the  chamber  is  proper 
mathematically  in  this  sense  is  likely,  though  not  com- 
pletely proV^tTi'  - 


t  ■'.  1. 1 ;,  i 

jHuu'iijijii  111 


-16- 

contour  la  represented  by  a  smooth  and  analytic  curve> 
However,  the  method  cannot  be  expected  to  be  applicable 
If  the  nozzle  Is  non-analytic;  e.g..  If  It  Is  pieced 
together  of  analytic  sections.   (The  approximations 
given  so  far  probably  remain  valid  If  the  curvature  of 
the  nozzle  contour  remains  continuous). 

More  serious  objections  can  be  raised  on  physical 
grounds.   It  Is  doubtful  whether  for  nozzles  with  wider 
divergence  the  assumption  can  be  upheld  that  viscosity 
can  be  Ignored,  or  that  the  flow  Is  steady  and  remains 
attached  to  the  wall.   In  this  connection  it  is  interest- 
ing to  confront  the  result  of  our  calculations  with  the 
experimental  results  of  R»  P.  Eraser  [3]. 

Eraser  has  tested  the  flow  of  compressed  air  into 
the  atmosphere  through  a  variety  of  nozzles  each  having 
the  same  expansion  ratio  (  I.e.,  exit  area  to  throat  area), 
six  to  one,  but  with  various  half -angles  of  divergence. 
He  has  employed  chamber  pressures  ranging  from  80  to  13 
atmospheres.   In  addition  he  has  made  shado.vgraphs  of 
the  emerging  Jet  and  has  measured  the  pressure  at  the 
wall  near  the  exit.  For  5°  and  10°  half- angles  of  diver- 
gence, the  pressure  agrees  rather  well  with  that  derived 
from  the  hydraulic  theory,  but  for  15°  divergence  the 
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pressure  is  about  5  to  8  percent  less,  while  for  20°, 
22-1/2  ,  25  ,  or  30   divergence,  the  pressure  is  10  to 
15  percent  greater  than  it  would  be  according  to  the 
hydraulic  theory.   The  latter  result  contradicts  the 
expectation  that  the  pressure  at  the  wall  would  be 
reduced  when  the  divergence  increases. 

Praser's  explanation  of  the  increase  of  wall 
pressure  for  wider  divergences  is  that  "it  may  be  that 
the  frictional  loss  in  the  convergent  portion  is  the 
more  i-:riportant  fact.   This  causes,  in  the  pressure  at 
the  throat,  a  consequent  drop  to  below  the  theoretical 
critical  pressure  and,  therefore,  a  consequent  drop  in 
the  pressures  along  the  divergence.   The  rising  mouth 
pressure  with  increasing  divergence  must  then  be  attri- 
buted to  a  further  departure  from  the  adiabatic  expansion 
in  the  divergence  due  to  general  turbulence  rather  than 
wall  friction."^  "''There  is,  however,  another  point  which 
is  important  and  probably  decisive  for  the  interpretation 
of  Eraser's  results.   The  contours  of  the  tested  nozzles 


*  A  similar  opinion  is  expressed  in  [4].   "it  will  be  seen 
that  the  losses  are  small  and  sensibly  constant  for  nozzles 
with  divergence  less  than  10°,  while  they  are  doubled  when 
the  cone  angle  changes  to  15°.   Since  in  the  latter  case 
the  length  of  the  nozzle  is  less  and  there  is  likely  to  be 
more  separation  (and  so  less  surface  friction)  it  v/ould 
seem  that  the  effects  of  turbulence  become  more  pronounced 
as  the  divergence  of  the  nozzle  increases." 
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consist  of  a  circular  entry  section  while  the  divergent 
section  issuing  from  the  throat  is  straight,  so  that 
the  contour  suffers  a  sudden  change  in  direction  at  the 
throat*   Certainly  this  change  in  direction  causes  a 
disturbance  which  should  be  considerable  for  wider 
divergences.   It  is  hard  to  say  of  what  this  disturbance 
consists.   The  flov;  will  probably  shoot  out  of  the  throat 
as  out  of  an  orifice  but  later  re-attach  to  the  nozzle 
wall.   In  any  case  a  shock  front  will  result  from  this 
disturbance  as  is  clearly  visible  in  most  of  the  photo- 
graphs of  the  emerging  jet.   Since  the  gas  suffers  a 
sudden  increase  in  pressure  when  it  crosses  the  shock 
front, it  is  possible  that  this  shock  is  mainly  responsible 
for  the  major  portion  of  the  10  to  15  percent  increase  in 
pressxire  observed* 

4.  Remarks  on  Jet  Detachment 

Jet  detachment  within  the  nozzle  has  frequently  been 
observed  to  occur  when  the  chamber  pressure  p   Is  so 
low  that  the  exit  pressure  that  would  result  from  undisturbed 
expansion  would  be  noticeably  less  than  the  outside  pressure 
p^  (here  assumed  to  be  1  atm).   In  breaking  away  from  the 
wall  the  jet  will  in  general  change  its  direction,  this 
change  being  accomplished  by  an  oblique  shock  front  beginning 
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at   the  point   of  detachment.      See  Fig.    2.      (One   might  expect 

vua.ll 
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Pig.    2 

a  state  of  rest  beyond  the  detached  jet,  but  Stanton^s 
experiments  indicate  back  flow  and  pressure  rise 
there).  Under  what  circximstances  jet  detachment  occurs 
depends,  as  A.  Busemania([  5]  pp.  395,399)  explains,  upon 
the  stability  of  the  resulting  flow  pattern  which  is 
essentially  determined  by  the  boundary  layer  at  the  wall. 
This  stability  may  be  very  feeble  as  certain  experiments 
of  A.  Stodola  ([6]  end  of  section  33)  indicate.   Strong 
dependence  on  Reynolds  number  is  indicated  by  experiments 
of  Stanton  ([1]  p.  324  Fig.  8)  who  has  tested  (  for  low 
pressure  ratios  p_/p„)  three  geometrically  similar 
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nozzles  and  also  varied  the  chamber  tempera ttire*  These 
and  other  experiments  indicate  that  the  point  of  detach- 
ment recedes  toward  the  throat  when  the  chamber  pressure 
p  (or  the  ratio  p  /p.)   is  lowered.  Praser  finds  in  his 

C  C   cL 

experiments  with  nozzles  of  various  divergence  that  jet 
detachment  Just  begins  when  the  exit  pressure  has  been 
reduced  (by  reducing  the  chamber  pressure)  to  about  .5 
to  .SS  atmospheres.   It  is  remarkable  that  this  should 
hold  rather  independently  of  the  length  of  the  nozzle. 
One  might  have  expected  that  for  a  widely  divergent  short 
nozzle  the  jet  would  have  broken  away  from  the  wall  with- 
out change  in  direction  as  soon  as  it  reached  outside 
pressure. 

In  any  case,  it  appears  that  the  jet  will  not  de- 
tach from  the  wall  before  it  has  reached  outside  pressiire. 
To  avoid  jet  detachment  one  should,  therefore,  make  the 
nozzle  so  short  tb^t  the  pressure  at  the  nozzle  wall  re- 
mains above  outside  pressure;  this  is  also  of  advantage 
for  producing  a  larger  thrust,  (see  the  following  section). 
For  a  widely  diverging  nozzle  the  hydraulic  theory  will  in 
certain  cases  predict  a  value  of  the  exit  pressure  above 
atmospheric  pressure  while  our  refined  formulas  indicate 
expansion  below  atmospheric  pressure  and  thus  the  possibility 
of  jet  detachment. 
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5.  Thrust 

For  the  reactive  thrust  produced  by  the  flow 
throxigh  the  nozzle,  one  can  easily  derive  a  simple 
approximate  formula  which  is  more  exact  than  that 
given  by  the  hydraulic  theory.   It  is  customary  to 
define  the  total  thrust  P  as  the  difference 

(5,01)  F  =  Fj^  -  P^ 

of  the  internal  thrust  F^^,  resulting  from  the  pressure 
acting  against  the  wall  of  the  chamber  and  the  nozzle, 
and  the  external  counter -thrust  F^  that  would  result 
if  atmospheric  pressure  were  acting  against  the  outer 
surface  of  the  body  in  which  the  nozzle  is  imbedded. 
(The  contributions  to  the  total  force  exerted  against 
the  body  due  to  the  deviation  of  the  external  pressure 
from  atmospheric  press\ire  are  considered  part  of  the 
drag  and  are  left  to  a  separate  investigation). 

To  evaluate  the  thrust  we  may  consider  the 
potential  surface  S^:  ^  =  ^^  =  const,  through  the 
exit  rim  of  the  nozzle  (Cf.  Fig.  3)  and  observe  that 
the  internal  thrust  P^^  (counted  positive  when  acting 
in  the  negative  direction,  i.e.,  against  the  stream)  is 
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equal  to  the  svm.   of  the  axial  component  of  the  monenttum  M^ 


I  aXm. 


I  atm. 


Fig.  3 


transported  through  S   to  the  outside  per  ixnit  time 
and  of  the  resultant  pressure  force,  P  ,  exerted 


against  the  surface  S   from  the  inside* 


(^"-) 


Since     vy/     is    the  mass  per  unit   time  flowing 

through  the  stream  tube    yr  ^  const,   and    try       is   the 

area  of  the  projection  perpendicular  to  the  axis  of 

that  section  of  S^  which  is  cut  out  by  this  stream 

e  "^ 


«•  I.e.  the  total  pressure  force,  P,-P  ,  exerted  against 
the  voltime  of  gas  inclosed  by  chamber,  nozzle  and  surface 
S  ,  equals  the  momentum  transported  per  unit  time,  M  • 


tube,  we  have 


M     =  «  /  q  cose  dy' 


e 


and 


P«  =  ■«  /  Pdy  . 


e 


6 


Thus  we  arrive  at  the  expression 


V  7 

(5.02)  F.   =  2tt  /  q  cobB  y/  df    +  2tt  /  pydy 

o  o 


the  integration  being  extended  over  the  surface  S  • 

6 


by 

(5.03) 


The  external  counter -thrust  is  obviously  given 


^a  -  Pa^ 


p   being  the  atmospheric  pressure  and 


(5.04) 


A  =  Try" 


being  the  area  of  the  projection  of  the  sxirface  S  on 
a  plane  perpendicular  to  the  axis.  It  is  convenient  to 
introduce   the   dimensionless    thrust  ratios 

(5.05)        K^  =  F^  /A,,J),,  ,    K^   =  F^  /A.„j),„.  ,    K   =  K^   -  K^ 
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where     A„      is   defined  by 


(5.06)  A.„.  =  tr  -V}^  =  iry^/  p,^,^     . 


The   significance   of     A^,.     is    that   it   is    the   throat 
cross-sectional  area  in  hydraulic  approximation 
(Cf.   Sec.    1). 

Introducing  the  variable  ^i    (Cf.(2.04)  )  along 
the  surface   S^  one  can  easily  derive  (Of.  Sec.  12) 
the  expressions 

(5.07)  K,  =  % «  /  (q/q^.  +  q«/q)  cose  Tid^ 

(5.08)  K  =  (P^/P,,)  (A/A„)  =  (pVp.J  {j/'Tlf     . 


To  evaluate  these  quantities  one  may  expand 
them  V7ith  respect  to  powers  of  7}   and  retain  terms  up 
to  the  third  order.   Instead  of  using  the  resulting 
expressions,  it  seemed  preferable  to  use  different  ones 
which  are  much  sin^ler  expressions  and  which  coincide 
with  the  exact  ones  up  to  terms  of  the  third  order  in  yi 
(Cf.  Sec.  12).   To  this  end  one  introduces  the  averages 
q  and   h   of  q  and  h   over  the  s\irface  S  S 

(5.09)    q  =  1  (qo,+  ^<l)  =  %tl  -^1  VS^^l 
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and 
(5.10) 


r^  =^  (h^  +  h^)  =h2[i+^(M2  .Dh^h;;'^'] 


derived  from  (2.13)  and  (2.15).  Here  q  and  h   refer 
to  the  nozzle  exit  rim.  Further  one  needs  the  expression 
for  the  half -angle  of  opening  at  the  nozzle  rim.  This 
is  given,  in  first  order,  by: 

e  =  h»77 

as  is  derived  from  (2.11)  and  (2.12). 

Employing  these  average  quantities  one  has  (up  to 
terms  of  third  order  in  Tj  ) 

(5.11)  K^  =  i-2  (q/q^.  +  q.^/q)  cos^  (6/2^ 

1  -^ 

(5.12)  K  =  (pVP-::-)  ii^  cos^  ^/2) 


(5.15)  K  =  [    ^  2  ^^/^-  •"  ^-^'^^  "  (Pa/P,Jh2]cos2(e/2) 
1  -/x 


This  approximate  expression  for  the  thrust  coefficient  is 
similar  in  form  to  that  of  the  hydraulic  theory, 

K  =  ^  i%M^^   +  %/%^   -   ^Pa/P*^^o  ' 

1  -^ 
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and  is  moreover  exact  for  pvirely  radial  flow  through. 
a  conical  nozzle. ^  ' 

The  hydraulic  expression  for  K  would  be 
approached  hy  making  the  nozzle  infinitely  long;  the 


<•   To  see  this  let  the  potential  surface  S   be  the 

^  e 

sphere  r  =  const.   Then,  in  view  of  q  =  q   and 


h  =  h   on   S  ,  one  has 


'n^  =   2r^(l  -  cos  e)h  '^ 


on   S^  and 
e 


/yj"2  /\os  e  dT?^  =  (1  -  cos  e)""^  /  cos  e  sin  ede=  cos^/2)^ 


hence,  from  (3.07), 


Ki  =r-H  <%/%  -^  ^^J%^   cosV2). 


1  - 


/^ 


Further,  from  y  =  r  sin  0,  one  has 


y^/07^=  h^  sin^  e/2   (1 


-  cos  e)   =  h^   cos^(e/2)^ 


whence 


K«  =  (P«/P*)  K  cos2(e/2) 


i\\ 
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maxlmvun  of  this  expression  is  easily  foixnd  to  be 

Therefore,  V^ax  '^'^^  serve  as  a  measure  for  the 
efficiency  of  the  exhaust  nozzle. 

In  the  design  of  a  rocket  exhaust  nozzle  it 
is  natural  to  aim  at  as  large  a  thrust  as  possible* 
The  problem  of  maximum  thrust  may  be  split  into  two 
parts: 

1.  An  indefinitely  long  nozzle  contour  is  given, 
which  may  be  cut  off  at  any  "end"  point.   The  thrust 
produced  by  the  remaining  section  is  to  be  investigated 
as  a  function  of  the  end  point.   The  maximum  thrust 

and  the  corresponding  "optimal" section  are  to  be  deter- 
mined. (Cf.  Malina  [7]). 

2.  A  "perfect"  nozzle  contour  is  to  be  found  whose 
maxlmuii  thrust  is  as  large  as  possible. 

•  Both  questions  can  be  answered  completely.   The 
answer  to  the  first  question  is  that  the  thrust  attains 
its  maximum  when  the  contour  is  cut  off  at  the  point  at 
which  the  gas  pressure  has  just  been  reduced  to  outside 
pressure.   In  that  case  the  thrust  is  given  solely  by 
the  moment\am  transport  per  unit  time,  (i.e.,  by  the  first 
term  In  the  expression  (5.02)).  For  the  hydraulic 

■■'4t 
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approxirtiation  this  fact  is  well  knovm  and  it  has  been 
frequently  assumed  tacitly.   It  should  be  emphasized, 
however,  that  it  holds  quite  generally.  '  This  state- 
ment can  easily  be  derived  from  relations  (3.07)  and 
'{3.08),  (Cf«  Sec.  12).   A  very  simple  proof  is  also 
given  in  an  earlier  memorandum  [  8] • 

As  to  the  second  question  a  rather  obvious  condi- 
tion was  derived  in  the  memorandum  quoted,  viz.,  that 
the  flow  direction  be  axial  on  a  sxirface  of  constant 
flow  speed  passing  through  the  rim  of  the  nozzle  mouth, 
provided  that  the  pressure  on  this  surface  equals  the 
outside  pressure.   How  to  construct  such  "perfect" 
nozzles  of  finite  length  will  be  discussed  in  the 
second  part  of  this  memorandum. 

The  task  of  constructing  a  "perfect"  exhaust 
nozzle  for  the  puj?pose  of  rocket  propulsion  is 
rather  academic,  since  the  maximtun  of  the  curve  repre- 
senting the  thrust  as  a  function  of  the  end  point 
is  very  flat  for  such  a  perfect  nozzle  contoiir.   The 
same  is  true  for  conventional  nozzle  contours o 

•5J-  P.  Malina,  [7],  in  an  attempt  to  take  the  divergence 
into  account  obtains  a  result  which  is  not  in  agreement  with 
the  above  statement.   The  reason  is  t-nat  no  account  is 
taken  of  the  changes  of  speed  and  press\ire  over  the 
cross-section,  although  those  changes  are  of  the  same 
order  of  magnitude  as  the  changes  in  direction  considered. 
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Experiments'  'with  more  than  a  hundred  tests  have  shovm 
that  the  thrust  Is  very  insensitive  to  the  length  of 
the  nozzle  and  that  Its  value  is  between  90  and  100  per- 
cent of  the  value  derived  from  the  hydraulic  theory. 
Hence,  shortening  conventional  nozzles  entails  "but  little 
loss  in  thrust  and  in  addition  keeps  the  pressure  well 
above  outside  pressure  thus  ensuring  the  absence  of  jet 
detachment  (Cf.  Sec.  4)» 

However,  the  thrust  is  more  sensitive  for  more 
widely  divergent  nozzles.   Therefore,  it  would  be  desirable 
to  obtain  information  about  the  tlirust  from  the  refined 
approximation  formula  (3. 13);  but  it  seems  difficult  to 
deduce  specific  properties  from  it.   Instead  we  shall 
formulate  certain  general  conclusions  that  were  gathered 
from  the  evaluation  of  the  thrust  in  various  examples 
discussed  in  the  following  section.   The  conclusions 
(of  necessity  somewhat  vaguely  formulated)  are: 

1.  The  maximum  thrust  of  a  contour  is  the  smaller 
the  more  divergent  the  contour  is* 

2*  The  maximum  thrust  obtained  from  the  short 
"optimal"  section  of  a  widely  divergent  contour  is  larger 
than  the  thrust  obtained  from  an  equally  short  section  of 
a  less  divergent  contour* 

"«•  B.  H.  Sage  [9] 
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The  latter  remark  Is  of  significance  since  it 
implies,  as  stated  In  the  introduction,  that  if  one 
wants  to  cut  off  a  section  of  the  nozzle,  and  Is  will- 
ing to  sacrifice  a  few  percent  of  thrust,  then  one 
could  obtain  the  same  tlirust  by  making  the  nozzle  still 
shorter  but  more  divergent*   It  should,  of  course,  be 
remembered  that  this  statement  should  be  qualified  by 
adding  that  it  holds  only  if  the  various  assvunptions 
enumerated  at  the  beginning  of  Section  1  remain  valid 
for  more  divergent  nozzles.   According  to  the  discussion 
in  Section  3,  it  is  doubtful  whether  this  is  so.   If 
the  assunrptions  are  valid  then  the  various  arguments 
discussed  are  all  in  favor  of  shorter  nozzles* 

6.  Examples 

As  a  first  example,  we  consider  a  set  of  nozzles 
which  are  approximately  symmetrical  with  respect  to  the 
throat.   These  contours  (referred  to  as  of  type  C-,)  are, 
to  some  approximation,  hyperbolas. (cf.  Table  I). 
They  are  obtained  by  setting 

(6.01)  1^0^^  ^  =  ^1  +  (  t/ef, 

e  being  an  appropriately  chosen  length.   The  function 


r  I  'i^ji  c^s  i^-^ii    iL.U 
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q  (  ^  )/q-;;.  is  then  to  be  found  by  solving 

-i-v 


1-/ 


~=  1  +  (|/e)2. 


Mq(  ^  )   is  then  determined.  With  eh^  =  (  ^  ^®y  A  "^^  ^/®)  » 

relations  (2.11)  and  (2.12)  becone 


e^h"  =  1 
o 


(6.02) 


7  =  '>?  /l  +(V/e)^  {  1  +  I  I^M^  -  1  -(^/e)2 


1  +(^/e)^^' 


To  first  approximation,  they  agree  with  the  relation 

hence  streamlines,  0^  =  const*,  and  potential  lines, 
X    =   const.,  are,  in  first  order,  confocal  hyperbolas 
and  ellipses  respectively.   In  any  case  the  streamlines 
approach  straight  lines  at  both  ends,  |  )j  |  -^  <'*=>• 

The  speed  q  along  potential  surfaces  \  =   const. 
is  given  by 

(6.03)    q  =  qo^  ^  ^  l^  -^  I  01/6)^(1  +  (|/e  )2|  . 
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Thus  the  relative  deviation  of  the  speed  from  the  value 

given  by  the  hydraulic  theory  approaches  zero  at  both 

ends  of  the  nozzle,  i.e.,  q/q^  -^  1  as  |  ^  |  -►  ©o  , 

A  set  of  streamlines,  potential  lines,  and  cxirves 

of  constant  speed,  for  /  =  1.2,  are  shovm  in  Table  I. 

In  Table  II  the  same  set  of  streamlines  is  shown  but 

so  magnified  that  they  all  intersect  the  throat  diameter 

in  the  same  point.   These  streamlines  then  represent  a 

set  of  nozzle  contours  with  the  same  throat  diameter. 

A  set  of  "curves  of  constant  thrust"  is  shov/n  in  this 

Table,  where  a  pressure  ratio  P-/p„  =  21.28  is  assumed. 

C    Si 

These  c\jj?ves  connect  the  end  points  of  nozzle  contour 
sections  for  v/hich  the  thrust  coefficient  as  calculated 
by  formula  (5.13)  is  the  same.   For  each  value  of  the 
thrust  coefficient  K  or  thrust  efficiency  K/K   , 
the  shortest  contoxxr  section  is  chosen  and  shown  separa- 
tely in  Table  IIA.   The  evidence  exhibited  in  Table  II 
clearly  confirms  the  statements  made  at  the  end  of 
Section  5. 

The  meaning  of  these  statements  will  become  clearer 
perhaps  when  contours  of  type  C,   are  compared  with  an 
essentially  different  type  of  contour.   As  such  we  choose 
for  a  second  example  the  ones  (referred  to  as  of  type  C«) 
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resultlng  from  the   assuriptlon 


(6.04) 


^o^n  =  1  +  ^\/^y 


c     "being   an  appropriately  chosen  length.      Here  we  find 
X  =    ^[i  -  [1  +    ii/cf]    (>l/c)^} 

y=    'yi[l-^    (^/c)2]     {l^^   [(m2    (^)    .l)(W(l)2).2(i)2J(^)^]. 

The  potential  lines  are  in  first  approximation  the 
circles  passing  through  the  points   x  =  +  c,7  =  0;  the 
streamlines  form,  in  first  approximation,  the  orthogonal 
set  of  circles.   The  speed  along  the  potential  lines  is 
given  by 

(6.06)   q  =  q^  (^  )  (l  +  [1  +  (^/c)^]  (^1/0)^]      . 

For  a  particular  case,  a  (nearly  hyperbolic)  con- 
tour C-^  and  a  (nearly  circular)  contour  Co  are  compared 
in  Table  III.   These  curves  are  chosen  such  that  they 

have  about  the  same  curvatures  at  the  throat.   More  specif i- 

2        2 
cally  we  have  chosen  e  =  6,  c  =  11  respectively.   The 

thrust  for  every  point  of  the  contour,  if  the  contour  were 

broken  off  there,  is  plotted. 


■^mejm 


f.^- 


iUyfi  i 


.u 


m%n^ 


-34- 

Two  observations  can  be  made.   Firstly,  the 
maximum  thrust  that  can  be  reached  at  all  Is  smaller 
for  the  more  curved  contour  Cg  than  for  the  more 
conventional  contour  C^.  Secondly,  the  values  of  the 
thrust  that  can  be  obtained  at  all  for  the  more  curved 
contou-r  Co  are  obtained  from  shorter  sections  than  for 
the  less  curved  contour  C-,* 

In  the  tv;o  examples  considered,  the  nozzle  con- 
tours were  approximately  symmetrical  with  respect  to 
the  throat  whereas  In  actual  nozzles,  the  entry  section 
is  usually  much  more  curved  than  the  exhaust  section. 
To  obtain  Information  on  flow  through  unsymnetrlcal  nozzles, 
contours  derived  from  (type  Cj) 


(6.07)     h^(  ^  )  =  ^(1  +  X^)  [!+(!+  >^f]    -X("l  +X) 

were  investigated.   The  slope  of  the  contoiirs  approach 
(approximately)  the  values  (v/ 1  +  X   -  X)''? 

and  -  ( \/l  +X   +  X)**?    as  ^  approaches  oo  and  -  oo 
respectively.   The  values  X  =  1.35  with  7l  +  X  =  1.67 
and  X  =  3.05  with  7l  +  X   =  3«20  were  chosen  for 
numerical  computation.  As  a  result  it  was  found  that 
the  flow  past  the  tliroat  Is  very  much  like  the  flow  past 
the  throat  for  an  approximately  symmetrical  nozzle  of 


-35- 


type  C, •  One  thus  gains  the  impression  that  the  entry 
section  of  the  nozzle  has  no  great  influence  on  the  flow 
past  the  throat»   It  was  also  found  that  the  convergence 
of  the  procedure  was  far  less  satisfactory  in  the  sub- 
sonic section  than  in  the  supersonic  section. 

Thus  the  examples  discussed  in  this  section  offer 
a  justification  for  the  general  conclusions  formulated 
at  the  end  of  Section  5. 
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PART   II 

ON  PERP^CT  EXHAUST  NOZZLES  AND  COMPRESSORS 

7.  On  Perfect  Nozzles, 

As  has  been  explained  earlier  (Sec.  5)  a 
nozzle  is  perfect  if  it  forces  the  flow  to  emerge 
from  the  nozzle  in  axial  direction  with  constant 
speed  and  a  pressure  just  equal  to  the  outside 
pressure.   The  reason  Is  that  then  the  reactive 
thrust  is  a  maxiniun>  i.e.,  greater  than  that  re- 
sulting from  any  other  nozzle  operating  with  the 
same  ratio  of  chamber  pressure  to  outside  pressure. 
This  statement  was  qualified;  it  is  valid  if  the  flow 
is  irrotational,  steady,  isentropic,  and  if  both  vis- 
cosity and  heat  conduction  can  be  ignored. 

Such  a  perfect  nozzle  can  be  designed  without 
difficulty.   As  a  matter  of  fact,  whenever  a  diverging 
exhaust  flow  is  given,  it  is  possible  by  re-routing  only 
a  certain  section  of  it,  to  make  the  flow  "perfect", 
i.e.,  to  guide  it  so  that  it  eventually  acquires  constant 
axial  velocity.  Every  streamline  of  such  a  perfect  flow 
yields  a  perfect  nozzle.   The  possibility  of  constructing 
a  perfect  flow  was  already  indicated  by  Prandtl  and 


mun^^ium 
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Busemann.  ' 


Perfect  nozzles  can  be  constructed  so  as  to 
produce  any  desired  exhaust  velocity;  that  is,  expressed 
in  dimensionless  terms,  the  Mach  number  of  the  exhaust 
flov/,or  what  is  equivalent,  the  ratio  of  chamber  to  ex- 
haust pressure  can  be  prescribed.   The  first  step  in 
the  construction  consists  in  securing  an  exhaust  flow, 

the  basic  "flow"  P  ,  which  leads  at  least  to  the  desired 

o 

exhaust  velocity.      To   this   end,    the  method  explained  in 

the   first  part  of  this  memorandum  may  be  employed.      The 

re-routing  process  for  a   two-dimensional  exhaust  flow 

Is   so   simple   that  it  may  be  described  briefly.     First 

the  point     A       on  the  axis   should  be  found  v;here   the 

basic  exhaust  flow     P^     attains   the  desired  exhaust 

o 


*  L.  Prandtl  and  A.  Busemann  [10]  p.  499;  Busemann  [5] 
p.  430,  and  [11]  p.  857.   In  the  first  two  of  these 
papers  a  graph  for  such  a  two-dimensional  flow  with  the 
Mach  number  M  =  1.85  is  shown  (Pig.  H   and  Pig.  51 
respectively).  A  similar  flow  picture  for  M  =  l.B  is  shown 
in  the  third  paper  (Pig.  10).   It  is  mentioned  there  that 
such  a  construction  is  also  possible  for  three-dimensional 
flow  with  rotational  symmetry;  a  procedure  is,  however,    •.. 
not  described. 

Nozzles  have  been  designed  on  the  basis  of  such  a 
construction  (following  the  two-dimensional  pattern)  by 
Mr.  Puckett  (Aero  Dept.,  California  Inst,  of  Tech.). 
Experiments  confirmed  that  these  nozzles  were  shock-free, 
(although  their  performance  otherwise  was  not  noticeably 
better  than  that  of  other  nozzles )[ 12]. 
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The  figure  shows  a  set  of  nozzle  contours  derived  by  magnifying  the  streamlines 

of  flow  C  .  End  points  of  sections  of  nozzle  contours  which  yield  the  sane 

thrust  are  connected  by  a  line.  The  theoretical  ratio,  kA   ,  of  this  thrust 

^  '  max 
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to  the  mximun.  thrust,  which  would  result  from  a   "perfect"  noz.le.   is   indicated. 
The  quantity    r^    here  equals  approximately  6  tir.es  the  curvature  of  the  contour 
at  the  throat  (cf.  p.  il). 
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Table   III 


The  figure  shov/s  two  nozzle  contours  of  the  type  C-,  and  Cp  so  chosen 

that  their  (approximate)  curvature  at  the  throat  is  the  same,  1/R  =  .20. 

The  thrust  ratio  K/K    as  a  function  of  a  variable  end  point  is  plotted. 

'  max 


U^ 


11  n«  &  nau-'^ri\ 


"fd*- 


«.- 


"wm 


« i^ 


ItU 


-38- 

velocity  q^.   Through  the  point  A   two  lines  are 

drawn:  the  baclcward  Mach  line  C   which  is  determined 

by  the  flow  P  ,  and  the  straight  line  D   which  vould 

be  the  forward  Mach  line  of  a  flow  E   with  constant 

o 

parallel  velocity  Qq  ,  (pressure  p   and  sound  speed 
c   being  determined  through  the  condition  of  isentropic 
change).   The  angle  between  the  line  D   and  the  axis 
is  just  the  Llach  angle  «^o   of  the  flow  E  .  Up  to  the 


Perfect  exhaust  flow 


line  C   the  original  flow  F   will  be  retained;  in 

the  sector  G  between  C   and  D   the  flow  will  be 

o        o 

changed,  and  beyond  D   it  will  be  parallel  with  the 

constant  velocity  q  •   To  determine  the  new  flow  in 

the  sector  G,   the  known  directions  of  the  Mach  lines 

of  the  flow  P^  should  be  marked  on  the  line  C^» 
o  o 
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Stralght  lines,  D,  should  then  be  drawn  from  C^  in  these 
directions  so  as  to  cover  the  sector  G.  Ptirther  the 
known  direction  of  the  flow  P^  on  the  line  C^  should 
be  marked  and  each  such  direction  is  to  be  transplanted 
parallel  to  itself  along  the  lines  D.   Thus  the  directions 
of  the  new  flow  P  in  the  sector  G  are  determined. 
Through  integration  of  this  field  of  directions,  beginning 

at  G  ,  the  streamlines  of  the  new  flow  F  are  obtained, 
o' 

Beyond  the  line  D   the  flow  is  to  be  continued  with 

constant  axial  velocity.  We  note  that  the  streamlines 

so  constructed  suffer  a  change  of  curvature  on  crossing 

the  Mach  lines  C   and  D  j  these  Mach  lines  will  there- 

o        o' 

fore  be  referred  to  as  "sonic"  lines. 

There  is  still  the  question  as  to  which  of  the 
streamlines  so  constructed  should  be  made  the  nozzle  con- 
tour.  A  reasonable  choice  can  be  based  on  the  following 
considerations : 

The  described  construction  is  possible  as  long 
as  the  set  of  Mach  lines  D  do  not  intersect  each  other 
in  the  sector  G.   This  might  happen  when  one  proceeds 
too  far  away  from  the  axis;  it  Is  certain  at  least  that 
the  angle  the  Mach  lines  D  make  with  the  axis  first 
increases  and  then  decreases  again  on  moving  away  from 
the  axis  along  the  Mach  line  C  •   The  place  where  this 
ang]HH|mmP^|VI|Bn  is  easily  determined;  it  is  the 
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Intersectlon  B   of  the  backward  Mach  line   C   with 
the  forward  Mach  line  which  begins  at  the  point  A, 
on  the  axis  at  which  sound  speed  Is  attained.   There- 
fore, one  Is  sure  that  the  construction  Is  possible, 
as  long  as  the  point  B  of  maximum  angle  is  not  passed. 
On  the  other  hand  one  wants  to  construct  the  nozzle  as 
short  as  possible.   Thus  one  is  led  naturally  to  choosing 
for  nozzle  contour  just  the  streamline  through  the  point 

The  construction  of  a  perfect  three-dimensional 

flow  and  a  perfect  three-dimensional  nozzle  (with  axial 

symaetry)  is  similar  to  that  for  the  two-dimensional 

flow  but  not  as  simple.   The  Mach  lines   D  are  no 

longer  straight  (except  D^)  and  the  flow  direction  is 

no  longer  constant  along  these  lines.   The  flow  is  to  be 

determined  from  the  basic  partial  differential  equations 

In  the  region  bounded  by  the  two  "sonic  fronts"   C   and 

o 

Dq»   These  equations  can  be  solved  by  a  numerical  procedure 
rather  easily  due  to  the  particular  type  of  Initial  condi- 
tion involved.   An  outline  of  the  procedure  is  described 
in  Section  13* 


%   This  is  apparently  the  choice  Prandtl  and  Busemann  [10] 
have  made  according  to  the  figures  they  show. 


-41- 

Buseraann  [11]  reports  about  perfect  flows 
determined  for  various  Mach  angles.   In  particular 
he  gives  values  for  the  raaxiniurn  half -angle  of  diver- 
gence along  the  nozzle  contour.   (This  half -angle  is 
the  angle  of  the  flow  direction  against  the  axis  at 
the  point  B   in  the  construction  above).   This 
maxiraum  half -angle  depends  on  the  Mach  number  M^ 
of  the  exhaust  flow.   It  becomes  zero  when  M  ->  1 

and  when  M^  -^  oo  ;  it  attains  a  maxiraum  somewhere 
o      ' 

between  M  =  3  and  M  =  5.   This  maximum  is  12°  for 
o  o 

the  two-dimensional  flow,  and  8°  to  9°  for  the  three- 
dimensional  flowf  (for  air  with  ^  =  1.405).   It  would 
also  be  interesting  to  knov/  what  the  length  of  these 
perfect  nozzles  would  be  for  given  throat  and  exit 
diameters.   Certainly  these  nozzles  would  be  rather 
long  and  very  softly  curved. 

Certain  limitations  for  the  curvature  of  the 
contour  at  the  exit  can  be  derived;  they  represent  the 
conditions  that  no  shock  develops  along  the  sonic  line 

D^«   Let  y  be  the  distance  of  the  exit  rim  or  lip 
o 

from  the  axis,  and  R  bo  the  radius  of  curvature  of  the 
contour  at  the  Up.   Then  for  two  and  three-dimensional 
flow 
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oc   being  the  Mach  angle,  connected  with  the  Mach  number 


)iFLf| 


M   of  the  final  flow  by  the  relation 
o 


I  =  1/3  in  «  • 
o    '      o 

For  the  derivation  see  footnote  in  Section  13,  p«  76  •  For 


M. 


=  3,  <?  =  1,4, 


we  find 


R/y  >  12.15. 


Also,  for 


M^  =  1.825 
o 


and  for 


R/y  >  5.7, 


M^  =  1.414 
o 


R/y  >  4, 8. 

The  significance  of  the  latter  value  is  that  for 
M  =  1.414  =  )J~2,    the  quantity  V      attains  its  mlni- 
num.   This  Mach  number  is  thus  most  favorable  as  regards 
the  lip  curvature  for  the  two-dimensional  flow. 

8.  Compressor  Flow. 

Isentroplc  flow  is  reversible.  Ifa  nozzle  is  so 
designed  that  a  shockless  exhaust  flow  through  it  is 
possible,  then,  in  principle,  the  reversed  flow  from 
exit  into  chamber  is  also  possible. 


/  f  P.  I  ,n  .I?  fl  ^  . 
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Thus  it  seems  that  all  that  has  been  said 
for  exhaust  nozzle  f lev?  would  hold  also  for  compressor 
flow.   It  would  thus  appear  to  be  possible  to  compress 
a  supersonic  stream  of  air  or  gas  so  that  it  comes  to 
rest  by  guiding  it  ttirough  a  reversed  exhaust  nozzle, 
a  "compressor".   In  the  following  discussion  we  present 
arguments  which,  with  certain  reservations,  speak  in 
favor  of  the  possibility  of  such  an  isentropic  compressor 
flow# 

Naturally  we  are  interested  only  in  isentropic 
compressor  flow  which  enters  the  compressor  with  constant 
velocity  in  axial  direction  and  hence  is  the  reverse  of 
a  very  special  type  of  exhaust  flow,  viz.,  of  what  was 
called  perfect  exhaust  flow.   Thus  the  compressor  should 
be  designed  as  the  reverse  of  a  perfect  exhaust  nozzle. 
V/hile  "perfection"  was  incidental  in  exhaust  nozzle  design, 
it  is  essential  in  compressor  design. 

The  exhaust  nozzle  is  designed  for  a  given  ratio 
of  chamber  pressure  p   to  exhaust  pressure  p  .   On  the 
other  hand  in  compressor  flow,  speed,  pressure,  and  density 
at  the  entrance  are  prescribed.  From  these  quantities, 
Bernoulli's  Law,  and  the  adiabatic  relation,  one  determines 
the  "stagnation"  pressure  p^.   corresponding  to  vanishing 


speed.  From  (1.03)  and  (1.05) 


r  T  ^  ^  -  1  „2  ,  >P  +  1 

The  stagnation  preastire  must  "be  maintained  In  the 
chamber  If  the  compressor  flow  Is  to  be  Isentroplc, I.e., 
P  *=?*••  Thus  the  ratio  of  stagnation  pressure  to 

C      So 

entrance  pressure  determines  the  design  of  the  compressor. 
The  question  remains,  however,  whether  the  flow  In- 
to a  perfect  compressor,  operated  with  the  appropriate 
chamber  pressure,  will  actually  follow  the  reverse  of  the 
perfect  exhaust  flow.  One  may  assert  that  this  is  so  if 
the  Influence  of  viscosity  and  heat  conduction  is  negligible. 
It  seems  well  established  that  the  Influence  of  viscosity 
on  Isentroplc  exhaust  flow  Is  very  smalli  '   One  should 
bear  in  mind,  however,  that  the  action  of  viscosity  is  not 
reversible.  While  for  the  exhaust  nozzle  the  rate  of 
growth  of  the  boundary  layer  developing  at  the  wall  is  less 
than  for  simple  parallel  flow,  this  rate  of  growth  is  greater 
for  the  compressor.  Since  the  flow  la  directed  against  a 
pressure  rise  there  is  the  danger  that  the  boundary  layer 
detaches  somewhere  in  the  diverging  section  between  throat 
and  chamber;  but  detachment  can  probably  be  avoided  by 
keeping  the  divergence  of  that  section  small  enough. 


«■  Except  that  viscosity  does  play  a  role  in  deciding  the 
position  of  shock  fronts,  if  they  occvir.  (Of.  Sec.  4,  p.  19. ) 
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Other  than  by  detachment,  the  boundary  layer 
will  hardly  influence  the  main  stream.   In  subsonic 
flow,  the  boundary  layer  influences  the  main  stream 
only  by  its  growth  inasnach  as  it  covers  a  gradually 
increasing  portion  of  the  main  stream;  this  portion 
will  be  very  small  if  Reynolds'  number  is  rather  large. 
In  supersonic  flow,  disturbances  in  the  boundary  layer 
are  propagated  immediately  into  the  main  stream  through 
the  Mach  lines.   It  does  not  appear  likely,  though,  that 
this  influence  of  the  boundary  layer  in  supersonic  flow 
is  strong  enough  to  cause  strong  changes  in  the  main 
stream.  Nevertheless,  this  question  needs  further  in- 
vestigation. 

Shock  fronts  are  another  disturbance  in  which 
viscosity  is  involved.   If  a  compressive  flow  cm^ves 
inward  too  strongly,  an  Isontropic  continuation  of  the 
flow  will  not  exist  and  a  shock  front  will  develop  in 
the  midst  of  the  stream.   The  shock  is,  however,  eliminated 
by  the  perfect  design  which  guarantees  the  existence  of 
an  isentropic  continuation  as  long  as  the  influence  of  the 
boundary  layer  does  not  interfere. 

In  expanding  flow  through  an  exhaust  nozzle,  strong 
shocks  can  be  created  by  sufficiently  increasing  the  pressure 
opposed  to  the  stream.   This  is  different  for  compressive 

SI 
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flow:   decrease  and  not  Increase  of  chamber  pressure 
will  cause  the  development  of  shock  front.  The  reason 
Is  that  the  stagnation  pressiire  decreases  when  the  entropy 
increases.  '  As  will  be  shown  in  the  following  section, 
it  is  necessary  to  operate  the  compressor  with  a  chamber 
pressure  less  than  the  "ideal"  one  for  which  the  nozzle 
was  designed,  since  otherwise  the  flow  would  not  be  stable. 

9.  The  Stability  of  the  Isentropic  Compressor  Flow. 

To  investigate  the  stability  of  the  Isentroplcally 
compressed  flow,  oiie  should  find  out  how  the  flow  would 
change  if  either  the  chamber  pressure  or  the  entrance  air 
speed  (but  not  the  compressor  design)  were  changed.   In- 
stead, we  ask  what  steady  flows  are  possible  under  such 
changed  conditions.   The  situation  Is  clear  for  an  Isen- 
tropic exhaust  flow.  liVhen  the  chamber  pressure  Is  changed. 


*  The  stagnation  pressure  is  given  by 

r  1  +'i     „2  ,N+  1   r    -5?i  ""* 
Pst  =  f  -rr  ^::-  ^        tP/°    ^ 

as  seen  from  (1.02)  and  (1.03).   The  critical  speed 
q^  is  known  to  be  unchanged  across  a  ehock,  and  pp 
depends  only  on  the  entropy  and  Increases  with  increase 
of  entropy. 
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the  flow  rate  will  change.   However,  when  the  outside 
pr3S3i..rv;  Is  changed,  the  flow  rate  will  not  change, 
but  the  jet  will  adjust  itself  to  the  varied  outside 
pressure  by  an  appropriate  pattern  of  shock  fronts 
and  rarefaction  waves. 

For  the  isentropic  compressive  flow,  the  response 
to  changes  in  the  entrance  or  chamber  conditions  is 
quite  different.   The  compressor  is  designed  to  be 
operated  with  a  certain  ratio,  p_/p-  ,  of  entrance 
pressure  p    to  chamber  pressure  p   and  a  certain 
Mach  number  M  ,  of  the  entrance  flow.   These  partic- 
ular  values  shall  be  referred  to  as  the  "ideal"  ones. 
Naturally,  we  keep  pressure  and  density  of  the  entering 
air  fixed.   The  quantities  to  be  varied  are  thus  the 
entrance  speed  q.   and  chamber  pressure  p  •   Apparently, 
one  obtains  the  clearest  pictxire  by  first  assuming  that 
the  entrance  speed  q_   is  changed  so  as  to  be  greater 
than  the  "ideal"  entrance  speed.   It  would  be  rather 
hopeless  to  appraise  the  influence  of  this  change  on  the 
flow  if  one  were  not  willing  to  consult  the  one -dimen- 
sional (hydraulic)  flow  theory  and  to  assxime  accordingly 
that  all  shock  fronts  occurring  in  the  compressor  are 
plane  and  cover  a  whole  cross-section.   Experience  with 
exhau^^MflaiHHlritLOw s  that  the  results  derived  from  such 


I 


Figure  R 


^■^ 


entrance 
Distribut 


tl:iroat 


chamber 


speed  and  pressure  in  a  coii5)ressor 
opLjmuxub  ^iT-xi  an  entrance  speed  greater  than  the 
ideal  one  and  with  various  chamber  pressures. 
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assvunptlons  represent  quite  well  the  cross-section 
average  of  wliat  occurs  in  reality.  "•'  The  hydraulic 
procedure^*''  then  indicates  that  the  flow  may  pass 
the  throat  in  a  supersonic  state  and  hence  be  a 
supersonic  expansive  flow  past  the  throat. (Cf •Pig* 5 
for  this  and  the  following  discussion),   jn  order  that 
the  flow  speed  reduces  to  zero  in  the  chamber,  a  shock 
front  may  occur  somewhere  between  throat  and  chamber 
(b,  c).  Where  this  shock  occurs  depends  on  the  chamber 
pressure.   When  the  chamber  pressure  is  increased  the 
shock  front  recedes  tov/ard  the  throat;  the  reason  is 
that  the  shock  then  becomes  weaker  since  the  entropy 
jvunp  decreases  with  increasing  stagnation  pressure 
(Cf.  footnote  on  p.  46  ).   The  chamber  pressure  for 


*  The  actual  pattern  of  shocks  is  more  complicated  and 
in  general  involves  jet  detachment.   That  the  sir^lifying 
assTimption  is  satisfactory  is  probably  due  to  the  re- 
attachment of  the  jet  -  a  fact  which  has  frequently  been 
observed.  Whether  or  not  such  re-attachment  will  occur 
in  compressive  flow  is  not  certain. 

•S-  As  seen  from  (1.08),  q./q..,.  increases  monotonically  with 

e   ■"""   p 

Mg  and  as  seen  from  (2.05),  h  (q^/q^  increases  monoton- 
ically with  q  /q^.j  hence  the  ratio  of  critical  to  entrance 
e  ■^>  p 

cross-section  area,  k,.Jk^  =   h  (q./q^, )  decreases  when 
%   =  ^e°e  increases.  When  q^  assumes  the  ideal  value  the 
tliroat  area  A.  equals  A„i  hence  A.>A3,  when  q   exceeds 
the  ideal  value.   In  that  case,  therefore,  the  flow  is  still 
supersonic  when  it  reaches  the  throat. 
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which  the  shock  just  stands  in  the  throat  (a)  may  be 

denoted  by  p^  =  p^  (q),   V/hat  happens  when  the 

chamber  pressure  Is  increased  beyond  p  seems  difficult 

c 

to  conjecture.   Probably  pulsations  in  the  flow  will 
result.   In  any  case  the  flow  will  be  unstable  when 

Pc  =  Pc- 

It  should  be  observed  that  there  is  another 

possibility  for  the  flow  when  p   <  p  ,  viz.,  a  flow 
v/hich  passes  through  a  shock  in  front  of  the  nozzle, 
(i.e.,  between  entrance  and  throat  (e)),  the  flow  then 
being  subsonic  when  passing  through  the  throat.   One 
could  imagine  that  such  a  forward  position  of  a  shock 
could  be  produced  from  a  flow  with  a  shock  placed  back- 
ward, by  sending  out,  from  the  chamber  against  the 
stream,  a  short  pressure  pulse  which  is  strong  enough 
to  push  the  backward  shock  over  into  forward  position. 
If  such  a  forward  shock  occurs  and  the  chamber  pressure 
is  decreased,  the  shock,  becoming  stronger,  will  move 
away  from  the  throat;  for  a  certain  chamber  pressure, 
it  v/ill  reach  a  limit  position  (f )  which  is  defined  by 
the  condition  that  the  flow  speed  is  just  sonic  at  the 
throat.  '  When  the  chamber  pressure  is  further  decreased, 

•K-  This  limit  position  advances  toward  the  entrance  when  the 
entrance  speed  is  increased.   It  is  interesting  to  note 
that  for  a  compressor  designed  for  the  Mach  number  Mg  =1.825 
this  limit  position  reaches  the  entrance  only  when  Mq  has 
been  increased  to  the  value  4.4,  as  is  easily  calculated. 
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the  shook  remains  In  the  limit  position  and  the  flow 

becomes  supersonic  past  the  throat  (g)#  The  section 

of  the  compressor  beyond  the  shock  front  then  acts  as 

a  forward  nozzle  and  the  adjustment  to  the  chamber 

pressure  will  be  achieved  by  a  second  shock.   It  is 

then  clear  that  to  every  chamber  pressure  below  p 

c 

there  are  just  two  flows,  one  being  supersonic  at  the 

throat  and  having  a  shock  in  backward  position,  the 

other  being  supersonic  in  front  of  the  throat  with  a 

shock  in  forward  position.   The  flow  with  a  shock  in 

forward  position  is  probably  the  more  stable  one.  At 

least  it  is  hard  to  conceive  how  a  shock  in  forward 

position  could  be  pulled  back  through  the  throat  by  a 

signal  sent  out  from  the  chamber. 

Suppose  now  that  the  entrance  velocity  is  again 

lowered  to  its  ideal  value.   Then  any  shock  in  forward 

position  or  at  the  throat  will  vanish  and  the  flow  will 

become  ideal.   Or,  if  the  chamber  pressure  p   is  less 
f  '  c 

than  the  ideal  one,  (also  denoted  by  p  ),  the  flow  will 
continue  supers onically  past  the  throat  and  adjust  itself 
to  a  chamber  pressure  less  than  the  ideal  one  by  a  shock 
in  backward  position. 

When,  further,  the  entrance  speed  q   is  decreased 
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below  the  ideal  value,  no  steady  flow  can  result,  what- 
ever the  chamber  pressiire  may  be.   The  reason  is  that 
the  critical  cross-sectional  area  A^,  as  determined 
from  the  hydraulic  theory  is  then  greater  than  the 
throat  cross -sectional  area  and  consequently  critical 
sound  speed  has  already  been  reached  before  the  throat; 
(Cf.  footnote  on  p.  48  )•   Insertion  of  a  shock  in  forward 
position  would  not  help  since,  due  to  the  increase  of 
entropy  across  the  shock,  the  critical  cross-section  would 
further  increase.  Probably  the  flow  will  become  oscilla- 
tory, and  if  the  chamber  pressure  is  high  enough,  the 
stream  might  even  reverse  its  direction  and  flow  out  of 
the  chamber  tlirough  the  entrance  opening. 

In  any  case, it  is  clear  that  the  ideal  flow  is 
unstable,  since  increase  of  the  chamber  pressure  above  p 
or  decrease  of  the  entrance  speed  will  cause  pulsations 
or  another  complete  change  of  the  flow. 

Therefore,  stable  operation  of  the  con^jressor  requires 
an  entrance  speed  q   greater  than  the  ideal  speed  for 
which  the  compressor  was  designed,  and  a  chamber  pressure 
p   less   than  the  maximum  value  P-Cq^)  which  corresponds 
to  the  entrance  speed  q  .   Thus  it  is  necessary  to  admit 
a  shock  in  the  operation  of  the  compressor,  but  the  shock 
strength  may  be  kept  small  and  its  position  can  be  confined 
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to  the  neighborhood  of  the  throat. 

10 •  The  Flow  at  the  Compressor  Entrance* 

Questions  of  particular  interest  arise  when  one 
considers  projectiles,  or  projectile-shaped  bodies,  which 
carry  a  compressor  of  the  type  discussed  in  Section  8. 
We  assume  that  the  compressor  is  bored  out  of  the  nose 
of  the  projectile,  as  indicated  in  Figure  6.  When  such 
a  projectile  is  opposed  to  a  supersonic  stream  of  air, 
the  air  has  the  opportunity  of  entering  the  compressor  and 


s\\ 


Fig.  6. 


of  arriving  in  the  chamber  isentropically  compressed. 
Whether  or  not  this  phenomenon  will  really  take  place 
is  the  question  to  be  discussed  in  the  present  section. 
(In  Section  11,  we  shall  discuss  the  question  of  the  drag 
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resulting  from  the  outer  surface  of  the  projectile). 

Naturally,  we  assume  tliat  the  coriipressor  Is 
a  perfect  one.  In  the  sense  of  Section  8,  since  -we 
know  that  otherwise  a  constant  supersonic  stream  could 
not  flow  without  interference  of  shocks.  But  even  if 
the  compressor  is  perfect  and  designed  for  the  Mach 
number  of  the  opposed  air  stream,  it  is  not  obvious  that 
the  air  will  enter  the  bore  with  constant  axial  velocity 
since  the  rim  may  cause  a  disturbance  of  the  flow. 

As  is  well  known,  ^when  a  projectile  with  a  coni- 
cal nose  moves  with  supersonic  speed,  a  shock  front  will 
develop,  which  is  either  of  conical  shape  and  begins  at 
the  tip  of  the  nose,  or  is  curved  and  stands  ahead  of 
the  projectile.   Which  of  the  two  cases  arises  depends  on 

the  nose  angle  and  on  the  Mach  nvunber  M   of  the  air 

o 

flow  relative   to  the  projectile. 

A  conical  shock  will  develop   if   the  nose  half -angle 
0      is   small  enough,    i.e.,    is    less    than  a   certain  extreme 
angle     ef^. ,   which  depends   on  the  Mach  number  M^.      For 
M_   =   1.91,    e.g.,    the   extreme    cone   angle    is      &?,    t.  =  39.5°, 

when  0     =  1.4;   when     M     -^  1,   6°   ^   -♦  0°,   and  when  M     •♦   co  , 

c  o 

0^    .    -»•    58  •      These  values    can  be  read  off  from  Busecoann's 
ext 

*  Of.    [13],    [14],    [15]    and  other   literature  quoted  there. 
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(D5]  Fig.  7  and  8)  or  Maccoll»s  diagrams  (Cf.[ 14] ,Pig.ll), 
Prom  these  diagrams  one  can  also  read  the  shock  angle  ^ 
that  corresponds  to  any  nose-angle  9^^  <  Q^^^   •  For 
M  =  1.91  and  0^  =  21°,  e.g.,  one  finds  ^^   =  40°. 


^% 


Fig.  7. 

In  case  the  nose  angle  Is  greater  than  the 
"extreme"  angle,  a  curved  shock  front  will  stand  ahead 
of  the  projectile.   The  position  and  the  shape  of  this 
shock  front  depends  not  only  on  the  Mach  nioraber  and  the 
nose  angle,  but  also  on  where  the  conical  part  of  the 
projectile  tapers  off j  and  quite  generally  on  the  shape 
of  the  projectile* 

When  the  projectile  carries  a  cocipressor  of  the 
type  considered  here,  the  exterior  flow  will  not  change 
very  much  except  at  the  rim.   The  flow  near  the  rim  is 
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locally  the   same  as   the    two-dimensional  flow  deflected 
by  a  wedge.      Accordingly,    the   condition  that  the   shock 
front  begins   at   that  rim  and  does   not   stand  ahead  of  it, 
is   the   same  as   for   the   two-dimensional  flow  with  reference 
to   a  wedge,   viz»,    that  the  wedge  angle     6      is   less   than 
a   certain  extreme  wedge   angle       6-_.4.   •      This   extreme  angle 


^% 


T^lg.  8. 

can  also  be  read  off  from  Busemann»s  diagrams.   It 
decreases  to  zero  when  M^,  decreases  to  1;  it  increases 
to  44,5  when  M^  increases  indefinitely;  its  value  for 


Fig.  9. 
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-se- 
ll =  1.91,    e.g.,    is  ©^=  21°.     Generally,    it   is    less 
than  the   extreme  angle  for  a   cone;    in  rough  approxiraa- 

*^°^     ®ext  =  -5  ®exf 

If  the  rim  angle  9^  is  greater  than  the  extreme 

angle,  the  shock  front  will  stand  either  ahead  of  the 
projectile  (Aa)  (Cf.  Figure  10)  or  ahead  of  the  rim  and 
be  drawn  into  the  projectile.   There,  either  a  disk- 
shaped  "Mach-shock"  perpendicular  to  the  axis ,  followed 
by  a  complicated  flow  pattern,  results  (Ab)  or  the 
shock  front  forms  a  double  cone  (Ac),  Which  of  these 
possibilities  occurs  depends  on  the  size  of  the  diameter 
of  the  compressor  opening  as  compared  v;ith  the  distance 
the  shock  front  stands  ahead  of  the  rim.   If  the  com- 
pressor opening  is  relatively  large,  the  drawn-in  shock 
front  may  fade  out  and  be  practically  sonic  even  before 
reaching  the  vertex  of  the  double  cone. 

An  entrance  shock  will,  however,  be  completely 
avoided  if  the  rim  angle  is  less  than  the  extreme  angle, 
(i.e.,  less  than  21°  for  11^=   1.91),  provided  the  entrance 
to  the  compressor  begins  like  a  cylinder  In  axial  direction, 
(  as  is  the  case  for  a  perfect  compressor);  then  only  a 
sonic  front  will  develop  from  the  rim  toward  the  interior. 
Thus  the  condition  is  found:  the  supersonic  stream  will 
enter  the  perfect  compressor  with  constant  velocity 
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Various    shock   conflgurat? ons  at   the    compressor  entrance. 
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wlthout  shock  Interference  If  the  rim  ar^^le  is  less 
than  the  extrerae  axigle  for  a  v/edgeo 

11«  Remarks  on  the  Drag  of  a  Projectile  Carrying  a  Compressor. 

Suppose  a  projectile  v/lth  a  conical  nose  Is  opposed 
to  a  supersonic  stream  of  air  so  strong  that  a  shock  front 
begins  at  the  tip  of  the  nose.   Then  the  pressure  on  the 
nose,  (whicl:i  is  constant  over  the  surface),  and  hence 
the  "pressure  drag",  i.e.,  the  total  force  resulting  from 
the  pressure  on  the  nose,  can  well  be  calculated  on  the 
basis  of  Maccoll's  oi'  Busemann's  diagrams.  As  Maccoll 
has  shown,  such  theoretical  results  agree  surprisingly 
well  with  experimental  results. 

If,  however,  a  compressor  is  bored  out  of  the  pro- 
jectile, the  tip  of  the  nose  is  replaced  by  a  rim  and  the 
flow  along  the  outer  surface  of  the  projectile  is  no  longer 
easily  determined.   Nevertheless,  some  rough  estimates  for 
the  pressure  drag  can  be  derived. 

It  was  explained  in  the  previous  section,  that  the 
shock  front  begins  at  the  rim  if  the  rim  angle  6   is  the 
extreme  angle  for  a  wedge,  Q^^^.,    (which  is  roughly  one- 
half  of  the  extreme  angle  for  a  cone,  Qly.^)*      Suppose 
that  this  condition  is  satlsfledj  then  a  shock  front  S 
which  behaves  like  that  formed  at  a  wedge  begins  at  the 
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w 


rlm«   The  angle  f3        of  this  shock  against  the  axis 
is  greater  than  the  angle  (8   of  a  shock  front,  S°, 
that  would  begin  at  the  tip  of  the  nose  of  an  unnruti- 
lated  cone  with  the  same  cone  angle  e»   (This  fact  can 
easily  be  read  off  frori  Busemann's  diagram).   Consequently, 
the  shock  S   at  the  rim  is  stronger  than  the  shock  S  • 
On  the  other  hand  one  should  expect  that  the  shock  front 
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Pig.  11. 

Plow  against  a  cone  with  a  bore. 
S  will  gradually  turn  over  into  the  direction  of  the 
shock  front  S°,  eventually  running  parallel  to  it» 

At  first  sight  one  might  also  expect  that,  along 
the  lateral  svirface  of  the  cone,  the  flow  will  approach 
the  flow  F   that  would  result  from  the  unniutilated 
cone.   This  assumption,  however,  would  not  be  correct. 
Since  the  shock  S  =  S^     at  the  rim  is  stronger  than  the 


uilUlhfK 


-59- 

shock  S^,  the  entropy  Increase  across  It  is  greater 
tlian  tlie  entropy  Increase  across  the  shock  S  ;  and 
both  entropy  increases  are  completely  independent  of 
the  cross-sectional  area  of  the  cylindrical  bore. 
Consequently,  as  can  be  shown,  the  sound  velocity  at 
the  lateral  surface  of  the  cone  is  larger  while  the 
air  velocity  there  is  smaller  tlian  the-  corresponding 
values  for  the  flov/  F°,  independent  of  the  size  of 
the  bore« 

If  the  conical  section  of  the  projectile  is 
sufficiently  long,  these  quantities  will  differ  from 
those  for  the  flow  F°  in  the  layer  whose  width 
approaches  zero*  As  in  the  case  of  a  boundary  layer 
induced  by  viscosity,  one  can  argue  that  the  pressure 
will  be  constant  across  such  a  layer*   Hence,  when 
the  size  of  the  hole  shrinks  to  zero,  the  pressure  will 
approach  that  for  the  flow  P°  while  the  values  of 
density  and  air  velocity  will  remain  smaller  by  a  finite 
amount  than  the  corresponding  values  for  flow  F  .   If 
the  press\ire  did  not  behave  in  the  indicated  way,  the 
"pressure  drag"  would  change  abruptly  when  the  hole 
disappeared;  this  would  be  absurd* 

Since  the  shock  S  s  S™  at  the  rim  is  stronger 

C  w 

than  the  shock  S   the  pressxire  p   Immediately 
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behind  It  is  greater  than  p°  behind  the  shock 
front  S°,  The  pressure  along  the  side  surface  of 
the  bored  cone  will  probably  decrease  mono tonic ally 
from  p  =  p  to  p  =  p°»  Generally  we  have  there 


w  ^      c 

P  >  P  >  P  • 


w       c 

The  values  of  p   and  p   can  easily  be  determined 

from  Busemann's  diagrams. 

The  preceding  remarks  lead  to  rough  estimates 
for  the  pressure  drag  against  a  section  of  the  cone  with 
a  bore*  For  a  complete  cone  the  pressure,  p  ,  Is  con- 
stant and  the  drag  is  simply 


d''  =  p""  A  , 


where  A  is  the  area  of  the  base  of  the  cone,  when 
the  cone  is  cut  perpendicular  to  the  axis  at  the  end 
of  the  section  considered* 
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Fig.  12. 
Flow  against  a  cone* 
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Consider  nov/  the   cone  with  a  bore  of  entrance 
cross-sectional  area      a. 


^\ 


Pig.  13. 
Plow  against  a  cone  with  a  bore 


w        c 
Prom  the  inequality  p  >  p  >  p   we  derive 


for  the  pressure  drag 


w 


p''  (A-a)  <  D  <  p   (A-a). 


If  the  cone  is  relatively  long,  ioe.,  if  a  <  A, 
the  value  of  D  is  probably  nearer  to  the  lov/er  limit. 
If,  hov/ever,  the  cross -section  of  the  bore  is  relatively 
large,  or  If  the  obstacle  deviates  soon  from  the  coni- 
cal shape  and  tapers  off  to  approach  a  cylinder,  then  the 
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upper  limit  p  (A-a)  Is  probably  a  better  approxima- 
tion to  the  coTT'ect  value  of  D, 


Pig.  14. 


As  an  example,  let  us  assume  that  the  incoming 
flow  has  the  Mach  number  1.91  and  that  the  cone  angle 
9  is  9.1  ,  Then  we  find  from  Busemann's  diagrams  that 
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Pig.  15. 
Plow  against  a  cone  with  a  bore 

the  shock  angle  for  a  wedge  with  6  =  9.1°  is  (9  "^  =  40°, 
while  the  shock  angle  for  a  cone  with  6  =  9.1   is 
(9°  =  32°,   Hence  the  shock  front  will  begin  at  the  rim 
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with  an  angle   of  40     and  eventually  turn  over   so  as   to 
make   an  angle  of  32     with  the  axis.      Further  one  derives 
from  Busemann's  diagrams 


pVPq  =  1.6     and     pVPq  =  1.19   , 


p  being  the  pressure  of  the  air  in  front  of  the  shock 
front. 

Hence  the  pressure  drag  should  lie  between 


1.2  p  (A-a)   and   1.6  p  (A-a) 


The  wide  margin  of  uncertainty  in  these  appraisals 
could  be  narrowed  down  only  by  a  more  refined  treatment, 
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APPENDIX 


12«  Mathematical  Details  to  Part  I  Concerning  Isen- 
troplc  Flow  Through  Nozzles* 

In  this  section  we  shall  supply  various  calcula- 
tions omitted  In  the  main  sections  of  Part  I. 

First  we  shall  give  a  detailed  derivation  of  the 
approximate  formulas  (2«11)  and  (2.15)  from  the  differ- 
ential equations  (2.01)  and  (2«02)» 

When  I  and  yj  ,   Instead  of  ^  and '\^  [Cf.(2.04) 
and  (2»05)],are  Introduced  as  new  variables  these  equations 
lead  to 


If  =  (qo/q)  cos  ®  » ||  =  -('n/y)  ^^  si^  ® 

||  =  (qo/q)  si-^  e  '  1^  =  iV/7)^^   cos  B, 


(12.01) 


as  Is  Immediately  verified*  The  Initial  conditions  for 
the  new  dependent  variables  are 

(12.02)      x=^   ,y  =  o   on  the  axis  77  =  0. 

Elimination  of  x  and  y  from  equations  (12.01)  yields 
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the  two  equations 

(    3(q^/q) 


(12.03) 


bnj 


3'>,  ^^^o^        at 


Si 


which,  together  with 
(12,03)„ 


|I   =  {'?7/7)h2  ^^g  Q 


serve  to  determine  q,  6,  y  from  the  initial  values 


(12.04)  q  =  q^{\)f   6   =  0,   J  =  0      on  the  axis»^  0. 


The  expansions  of  these  quantities  with  respect  to 
powers  of  07  have  the  form 


'    7  =  7^+7307^+...   J 


(12.05) 


e  =  e-,07  +.. 


1'^ 


q/q^  =1  +  80'^   +. 


Prom  the  last   expansion  one  derives,  by  (2.06)   and  (2.07) 
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1     +    7T     [M 


-     1]    Xg^?^     +•••     I 


and  by   (1.05)    and   (l.08). 


p  =p^{l  -    J?m2   S2'>?^  +...} 


Insertion  in  (12.03)   yields 


r  71  =  h^( n 


(12.06){ 


j3=-5  t^o   ^^)-  I3h^(ni^j   (^)--5[i^i(^)    ]' 


For   the  expansion 


X  =  I  +  Xg  -yj      +, 


one  then  finds  from  bTir/dyi=   -  (  'V}/y)h  sin  6, 


Insertion  in  (12.05)  gives  the  expansions  (2.11)  to 
(2.15). 

Secondly,  we  shall  supply  the  steps  necessary  for 
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the   calculations   of   the  thrust.   Section  5.     We 
Introduce  the  variable    17    In  the   Integral   (5.02), 
which  is   to  be   extended  over  a   line     i    =   const., 
and  employ  the  following  relations* 


Yd>=  |O^q,,.07d')7 


from   (2.04),     and 


2  2 

ydy  =  h     cos     9  7]d7| 


=  f*.-.%.  /°"  q"      °02   e  r^drf 


from  (12»03)^  and   (2. 06). 


Furthermore,   we  use 


(12.07) 


-1  12        1   -  V      2 


(12.08) 


/o"     p   = 


1        r    2  ,2, 

1  +  >J       ^' 


from  (1.02),    (1.04)    and  c^  =  q.. 


and 
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Thus 


K     =  i-t^  -^    /  (q/q.)    cos  6  fi^yj 
■"•       1  -  V       *>?      o 


+ —-2   /  C  q^j/q  -  T»q/q^  1  cos  e  yid-rj 

1    -   V  ^  o 


12^ 


1  -'J      ->?"   o 


Using  the   expansions    (2.13)    and   (2«14)   up   to  second 
order   in    07   ,    one  obtains 


i  =  7-J7  Hq^q.^  +  q^-/qo)(i  -  h^l^  '»?^) 


K,    = 


+  ^qoA-  -q../qo)  VS^^  3  • 


It  is  easily  verified  that  the  expression  given 

in  (5© 11)  possesses  the  sane  expansion  up  to  second  order 

p 
in  77   (or  up  to  third  order,  since  no  third  order  terms 

occur). 

The   expansion  of  the    expression     K       given  by 

(5. 08),    is   fo\ind  from   (2*12), 
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It  agrees,  up  to  second  (or  third)  order,  with  the 
expansion  of  the  expression  given  in  (5. 12),   Thus 
the  basic  approximate  expression  (5.13)  for  the  thrust 
is  established. 

A  proof  may  be  given  that  the  exact  thrust  F  as 
a  function  of  the  cut-off  point  of  a  given  contour  is  a 
maxim\xn  when  the  pressure   at  the  cut-off  point  coincides 
with  the  outside  pressure.   To  this  end  we  make  use  of 
identities 


If-  =  2(qyq)  y  sin  0 


S| 


and 


aT  T-TT  ^^^'^"-  *  '^"■^'^^  cos  eo,  =2  ^(p/p^,)(q/q„)y  sin  0, 

which  can  be  verified  by  using  relations  (12.01),  (12.02), 

(12.07)  and  (12.08).   Integration  with  respect  to  t^  yields, 

using    (5.07)    and   (5.08) 

3^  _  2pq^y  sin  9 

di        i  P.;.q 


a  _  gp^q^  7   sine 
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where  it  Is  understood  that  the  Integral  Is  evaluated 
at  the  nozzle    contour* 


Hence 


,._Jl;M^^]  ,3.ne. 


d\  \    P 


« 


AS  a  result  it  follows  that  K  is  stationary  when 
at  the  cut-off  point  either  p-P^  or  9=0,  the  latter 
condition  would  necessitate  cutting  the  contour  at  the 
throat.  That  K  assumes  a  maximum  in  the  former  case 
and  a  minimum  at  the  throat  could  easily  he  shown.  We 
refer  to  a  different  presentation  in  an  earlier  memoran- 
dum [8]. 

13.  non»truction  of  ^  P^T-fect  Three -Dimensional  Flow  from 
Any  Expanding  Flow  with  A^lal  S^metry. 

It  was  shown  in  Section  7  how  any  expanding  two- 
dimensional  flow  with  axial  symmetry  can  be  turned  into 
a  constant  parallel  flow.   In  this  section  we  shall  explain 
how  the  same  problem  can  be  solved  for  three-dimensional 

flow. 

suppose  the  flow  is  to  emerge  with  the  speed  q^. 

Then  one  determines  the  point  A^  on  the  axis  at  which 
this  speed  is  reached.   Through  this  point,  the  Mach  cone 
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C   Is  drawn  In  upstream  direction.   The  original 
flow  is  then  to  be  retained  only  up  to  this  cone  C  • 
Further  one  draws  through  this  point  A   in  down- 
stream direction,  the  straight  cone  D   whose  angle 
against  the  axis  is  the  Mach  angle  belonging  to  the 

speed  q  •  Along  this  cone,  all  quantities  character- 
o 

izing  the  flow  are  constant;  their  values  are  there- 
fore the  same  as  at  the  point  A  •  Since  these 
quantities  are  given  along  both  cones  C   and  D  , 
the  flov/  in  between  can  be  determined  by  solving  the 
partial  differential  equations  of  the  flow  problem* 
Mathematically  speaking,  the  problem  is  an  "initial 
value  problem"  for  a  hyperbolic  differential  equation, 
the  initial  values  being  given  on  two  "characteristic 
lines. "^""^ 

To  this  end  the  basic  differential  equations 
(2.01)  and  (2.02)  will  be  re-written  by  referring  them 
to  characteristic  parameters^'''*' i.  e. ,  to  variables  which 
are  constant  along  the  Mach  lines.   Instead  of  the  speed 


*  Cf.  these  notions  [2]  Ch.  5  and  [16], 


^c-  The  method  explained  in  the  following  can  also  be 
employed  for  a  variety  of  other  three-dimensional  super- 
sonic flow  problems. 
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q     we   introduce   the  angles    7j ,     S ,   and    t    by 


(13.01) 


tan*);      =     0 


M^   -  1 


(13.02) 

and 

(13.03) 


tan  S      =      utanT) 


t,    =  ^*S  "  n^        • 


Differentiating  (13.03)  we  obtain 


.        p  p 

dc,  =  ucos  S  d  tanS  -  cos  TjAtano^ 


2        2 
(cos  S  -  cos  0^  )d  tan'y^ 


or  since  by  (1.08) 


d  tan  yi 


1 2d£ 

,T    2,    2f    2      q 
(1  -yu.  )cos  6  cos  nj  ^ 


(13.04) 


d^  =  tan^^  — S 

q 


and  by  (1.07) 


(13.05)   ^L£^ 


=   -  tan  T?  dt   • 


It  is  clear  that  q  may  be  considered  a  function 


^^iiLHuull  ill 
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of  4  ^^d  consequently,  also  yf ,    S   ,    o,  p,   and  /j    may 
be  considered  functions  of  ^  •   The  latter  three  quantities 
are  easily  expressed  in  terms  of  6  •   This  is  seen  from 
the  relation 


2f 
cos  6   = 


2    2  2 


which  by  (1.04),    (l.05)    and  (l.06)    leads   to 


2  V  +  2,  2  V 

(13.06)       £-  =  cos5     ,  ;2.  =  cos  S  ,  -£.  =  cos     5  , 


Purther,  one  easily  verifies 


(13,07) 


q  _  sin  i 


^3in>j 


cos  i 

COSIr^ 


When  we  now  eliminate  0     or  y'  from  equations 
(2.01)  and  (2«02)  by  expressing  the  identities^  ' 


a       -  rt   =0.  (Mr^)       -  i^   )^   =  0  In  terms  of  the  right- 
'"xy   "yx   ^'  ^  y  'xy  ^  r    'jx 


hand  members  of  these  equations,  we  obtain  a  simple  form 
of  the  result  by  virtue  of  (13.04)  and  (13.05) 

^   Subscripts  X  and  y  indicate  partial  differentiation 
with  respect  to  these  quantities. 


'    .,    f     :    f: 


ii^mmififh 


(13.08)' 


-74- 


[cot  '>?4^-©_]co9  e  -  [cot  y)(^+  ©^]sin  e  =  0 


[tan  rjd^-e^] sin  e   +  [tan   174    -6    ] cos   e-y'-'-siii  0=0 


X  r 


To  handle  o\ir  problem  It  I3  natural  to  introduce 
characteristic  parameters  c    and  t  since  then  the  two 
boundary  cones   C   and  D   are  given  by  fl"  =  const* 
and   T  =  const,  respectively.  We  assume  in  particular 
that  C   is  given  by  c  =  0  and  D   by  r  =  0. 


¥  X 


Fig.  16. 


Tlie  characteristic  directions  are  then  characterized 
by  the  condition  that  they  make  the  Mach  angle  oc  with 
the  flow  direction  or  the  angle  7^  v,/ith  the  normal 
to  that  direction.   We  count  O"  and  T  positive  away 
from  the  axis  as  indicated  in  Fig.  17. 


y 

\ 

\ 

y 

^^ 

T 

A 

^ 

-^\ 

\ 

Fig.  17. 

Instead  of  derivatives  k—  and  -^       we  employ 
differentials   d   and  d   for  the  corresponding  directions. 
The  equations  for  x  and  y   then  are 


(13.09) 


^s^   =  S  d^y,  d^x  =  -  Td^y 


with 


(13.10) 


S  =  tan  (0^  -  9),  T  =  tan( '>^  +  0). 


Equations    (13.08)   are,   by  virtue   of   (13.09),    equiva- 
lent with 

©)   =  Pd^y,   d^(^  +  9)   =  -  Qd^y 
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where 


/ 


P  = 


(1.3.12)/ 


_1  cos  7}  sin  Q 
^  cos(  '»7  -  e) 


-1 


1 

y 


cot  e  +  tanT] 


=  y   sin  e  [cos  6  -  S  sin  6] 


Q  = 


1   cos  7}  Sin  e  _  _1  1 

•^  cos(')j  +6)    ^   cot  e  -tan-y^ 


\ 


-1 

=  y    sin  e  [cos  e  +  T  sin  9] 


Relations  (13.11)  can  easily  be  verified;  their  form 
exhibits  the  fact  that  (l3»09)  determines  characteristic 


directions 


(-) 


*  From  the  formulas  (13.01)  to  (13.12)  one  easily  derives 
the  upper  limit  for  the  entrance  curvature  as  given  in  Sec. 7. 
The  relation  which  one  obtains  from  (13.09)  by  elimination 
of  X  reduces  to 

(1)     2  tan  rjod^^   y  +  sec^  T7od^(  tj  -  e)d^y  =  0 

on   the    line     Dq*      From  d^   (cos^S      -   cos^ ')7  )d  tami  = 

(1   -  yu^)    cos' 

and  d^^  =  -  d^e  on  Dq,  one  derives  d^(e 


•^  h     tan^i;  dT^  or  d^^  =  (1  -  /i?)"-^(cot^'>j  + /x^;d^. 


-7,)  ={1-J\c 


JSC^TJq 


where  the  parameter  T  is  so  chosen  that  d^9  =  1  on  Dq.  If 
the  Mach  lines  D   converge  at  the  point  Aq  we  have  d^y  =  0 
there,  and  relation  (1)  can  be  Integrated  to  give 
(2)  2  tan  T7odT.y  =  (1  - /^2)-l(jgj.2  ^  ^gg^2  ^^  y  _  f-Q  y. 

On  the  other  hand,  along  a  streamline  through  a  point  of 
Dq  we  have  dy  =  0,  hence  ds  =  -dx  =  2  tanTJod^yj  further 
d  e  =  d^G.   Therefore,  the  radius  of  curvature  is  R  =  ds/d  9 
2  tan  -y^Qd^y  =  Ppy.   If  the  Mach  lines  do  not  converge  on 


D. 


,^^bvlous'j 


ive 


R  >  To  y. 


ii 
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From  relations  (2#01)  and  (2»02)  together  with 
(13.06)  and  (13.07)  one  obtains 


(13.13)    d^^  = 


%^       s  in  8 


^  cosC^  -  e) 


,  ,     ,  %'<r        sinS 


^    003(1^  +  e) 


■1^7 


(13.14)  d,/.  ^^,^,  =o£llii_4^y2,^^^2  ^  ^^^^^  S2S^^!-^  a^y^ 


cos('>j  -  e) 


cos(n  +  e) 


As  initial  conditions  we  have 


(13.15)  y  =  cot  ti^ix   -  x^),  e  =  0,  i  =  4^,  ^  =-£ 


q.,  sin  S 


•K-  COS   01 

'  O 


2        cos^^-^l  g 


COS 


'^o 


-  7  » 


where  "Vj  ,    S^,  and   <^   are  determined  fron  (13.01), 
(13,02),  and  (13.03)  for  q  =  q^. 

Along  C   we  have  the  condition  that  ^,  y,  6,  q 
and  hence  y\ ,     8,4  ^^®  those  functions  of  x  and  y  that 
characterize  the  incoming  stream.  A  further  relation 
between  x  and  y  (  if    and  9)  is  given  by  the  second  rela- 
tion (13.09)  . 

V/e  note  that  the  differential  equations  (13.11) 
are  rather  singular  at  A  .  On  D  ,  i.e.,  T  =  0,  we  have 


A 


V 


J^' 


ii  1  1^  3    /r  i 

mi 
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e  =  0  and  hence  p  =  Q  =  0.   On  C^,  i.e.,  (T  =  0,  the 
ratio  ^   ■   approaches  a  value  different  from  zero, 
hence  P  and  Q  do  not  approach  zero  as  T  ♦  0  on  ff  =  0. 
Nevertheless  it  is  likely  that  the  problem  has  a  unique 
solution. 

To  obtain  the  solution  numerically  we  intend  to 
employ  the  method  of  finite  differences.  We  place  a  net 
of  lines  0"  =  const,  and  T  =  const,  across  the  sector 
between  C   and  D^.  By  A,  and  A_  we  then  indicate  the 
differences  of  values  at  neighboring  net  points  in  the 
<r  and  t  direction,  respectively.  More  specifically 
let  z  be  a  quantity  at  one  net  point;  then  we  denote  by 
z*  and  z"  the  values  of  that  quantity  at  the  next  net 
points  reached  by  going  backwards  in  ff-  or  t  direction 
respectively. 


z* 


->-  (T 


Pig.  18. 


-■■:   S  I 
If  \ 
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We  then  set   A^z  =  z  -  z"*"   and   A^  =  z  -  z' 


( 13, 16 ) 


1 


S  = 


I  (S"^  +  S),  T 


I  (T"  +  T) 


P  = 


I  (P"^  +  P),  Q  =  -|  (Q"  +  Q) 


Then  we  replace  the  differential  equations  by 
the  difference  equations: 


(13.17) 


A^x  =  S  A^y  ,   A^  =  -  T  A^y 


(13.18) 


A(4  -  e)  =  P  A^y  ;  A^(4  +  9)  =  -  Q  A^y 


(13.19)  A0=^f  "^^^ 


•^   -^  1003(7?  -  e)/ 


"  \cos(')7  +  e)/ 


(13.20)  A;v^=  /«„q.:. 


cos      b 


COS  (07  -  e), 


A   2   .^2   ^    /cos^"^  •*••'•&  \.   2.   . 
A^y  ;  a;^     =  f^.(LJ 7- —  A^y 

1003(07  +  e)/ 


We  eliminate  x  '  from  equation  (13.17)  and  obtain 


(13.21) 


y  =  (x  -  X*  +  S  y*  +  Ty  )  /(S  +  T)  . 


Further,  from  (13.18)  we  have 


(13.22)    24  =  4*  +  <  "  -  (e"*"  -  ©')  +  P  A,y  -  Q,£^^  y 


i.. 
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(13.23)   2e=e'^  +  e"-(4'*"-  4")  -PA^y-QA^y 

The  procedure  is  then  as  follows:   Suppose  the 
values  for  (+)  or  (-)  are  known;  then  one  first  replaces 
5,  T,  P,  Q,   by  S**",  T",  P"*",  q"  respectively  and  determines 
y  from  (13.21),  x   from  (13.17)  or  (13.18),  t,    and  9 
from  (13.22)  and  (13.23).   Prom   ^  one  is  to  find  'v^   • 
To  this  end  one  may  employ  a  diagram  or  a  table  by  Busemann 
([5],  pp.  424  and  425)  provided  one  is  v»illing  to  set 
%   =   1.405.   How  to  do  this  may  be  shown  by  an  example. 
Suppose  4  =  21°  is  found.   Then  one  follows  the  epicycloid 
in  Buseraann's  diagram  [5]  until  Its  polar  angle  is  21  • 
The  corresponding  point  has  a  distance  from  the  origin 
which  belongs  to  the  "pressure  number"  979.   In  Buseraann's 


Fig.  19. 
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Table  I  one  then  reads  off  q/c  =  w/q  =  1.815.   This 

ntunber  Is   sec  nj     hence  07  can  "be  determined*   Having 

foiind  97  one  can  determine  S,  T,  P,  Q«   Now  one  is 

able  to  calculate  the  average  values j  and,  employing 

them,  one  can  repeat  the  process  until  the  results  do 

not  change  appreciably.   The  scheme  corresponds  to 

Moulton's  method  for  solving  differential  equations.   It 

seems  necessary  to  follov;  this  complicated  scheme  in  the 

present  case, at  least  near  the  point  A  ,  In  view  of  the 

strong  singularity  there. 

Since  all  desired  quantities  are  given  on  C^  and 

D   it  is  clear  that  the  solution  can  be  determined  in 
o 

all  net  points  within  the  sector  between  C   and  D  • 

Example 

As  an  example,  we  consider  the  case  that  the  incoming 
flow  is  purely  radial.  Such  a  flow  is  cliaracterized  by 
the  solutions 

(13«24)  X  =  r  cos  8,  y  =  rsinej 

where  r  =  r(q)  is  given  by 

(13.25)        r  =  r^,  ^  q^  Pit /^C       =  ^^j>> 

r„   is  the  distance  from  the  origin  at  which  the  critical 


K'jvij^g!, 


IhVI  n 
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speed  Is  reached;   q/q^.  ,  p/  f>y     are  given  by  (15.06) 
and  (13.07)  •   One  easily  verifies  that  relations  (13.08) 
are  then  satisfied,  if  one  observes  that  r"  dr  = -^  tan  Vj  d^ 
holds  by  virtue  of  (13.25)  and  (13.05).  Also  we  find 


(13.26) 


-yr^   =  2r|^^q^,(l  -  cos  e  ) 


and 
(13.27) 


ij  =   /  q(r)dr  • 


To  obtain  the  equation  of  the  characteristic  D   we 
insert  (13.24)  in  the  second  equation  f 13.09)  and  find 
r"  dr  =  -  tanTj  d  0.  Comparison  with  r"  dr  =  w  tan 'xj  d4 
yields 


(13.28) 


e  =  i(  4„-<:). 


(^   being  the  value  of  ^  at  the  point  at  which  the 

characteristic  intersects  the  axis,  (jt  is  interesting 

to  note  tliat  the  equation  of  the  characteristic  for 

two-dimensional  flow  v/ould  be  6  =  4^-  4  i  for  three- 

o    ' 

dimensional  flow,  therefore,  the  same  state  as  for  tv/o- 
dlmenslonal  flow  is  reached  on  a  characteristic  on  turning 
tl23Raii^^3^1|j^he  angle.  )  To  a  given  value  of  ^  one 


mm 
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determines  4  >  then  reads  off  from  Busemann's  diagram 
[5]  the  "pressure  nvanber"  and  takes  the  values  of 
^%   iO.,j./q|0    which  give  r/r..,,  by  (13.S5).  All  other 
quantities  are  then  readily  obtained* 


\  \    \  \ 


.<V^**^^  .   .       ^5-H     ^.    s^'V.., 
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can  be  used.  Braun' s  method  is  carried  through  for  the 
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